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ABSTRACT 


Measurements  of  the  boundary-layer  characteristics  and  pressure  distribu- 
tions on  three  axisymmetric  bodies  with  and  without  a propeller  in  operation 
are  correlated  with  analytical  techniques  for  computing  ( I ) the  interaction 
between  the  hull  boundary  layer  How  and  the  potential  flow  and  (2)  the 
interaction  between  the  operating  propeller  and  the  thick  stern  boundary 
layer.  The  agreement  between  the  measured  and  computed  shear  stresses, 
pressure  distributions,  and  velocity  profiles  was  satisfactory  except  over  the 
last  5 percent  of  body  length  on  the  two  fullest  afterbodies.  An  inviscid 
propeller-stern-boundary-layer  interaction  model  was  developed  and  found  to 
give  accurate  predictions  of  effective  velocity  profiles.  The  experimental 
results  show  that  the  potential-flow  propeller-hull  interaction  methods  give 
accurate  predictions  for  thrust  deduction  and  the  propeller-induced  pressure 
distribution.  The  present  analytical  techniques  will  be  useful  for  predicting 
the  full-scale  effective  velocity  profile  for  propeller  design  when  account  is 
taken  of  the  effects  of  Reynolds  number,  roughness,  and  propeller  suction  in 
modifying  the  nominal  velocity  profiles. 

ADMINISTRATIVE  INFORMATION 

The  initial  experimental  work  reported  herein  was  supported  by  the  independent 
exploratory  development  program  ( 1 1£ L> > David  W.  Taylor  Naval  Ship  R&D  Center 
(DTNSRDC).  Funding  was  provided  under  Element  62766N,  Work  Unit  1-1552-119.  All 
subsequent  experimental  and  analytical  investigations  were  supported  by  the  Naval  Material 
Command  direct  laboratory  funding  under  Element  62543N,  Work  Unit  1-1520-004. 


INTRODUCTION 


The  radius  of  many  ship  propellers  is  about  the  same  order  of  magnitude  as  the  thickness 
of  the  boundary  layer  at  the  propeller  location.  In  these  cases,  the  (low  in  the  vicinity  of  the 
propeller  is  characterized  by  the  presence  of  a thick  and  possibly  separated  turbulent 
boundary  layer.  An  operating  propeller  produces  an  upstream  suction  which  results  in  an 
augmentation  of  stern  pressure  drag  and  skin-friction  drag  due  to  the  acceleration  of  (low. 
Naval  architects  refer  to  this  drag  increase  due  to  propeller  suction  on  the  afterbody  as  thrust 
deduction;  it  is  usually  determined  by  conducting  conventional  resistance  and  self-propulsion 
model  experiments  in  a towing  tank. 

The  suction  of  the  propeller  also  increases  the  velocity  in  the  entire  boundary-layer 
region.  The  effective  inflow  velocity  experienced  by  the  propeller  is  the  result  of  the 


interaction  of  the  propeller  and  the  thick  stern  boundary  layer,  and  it  is  different  from  the 
nominal  velocity  distribution  measured  in  the  absence  of  the  propeller.  The  nominal  velocity 
distribution  at  the  propeller  plane  of  the  ship  model  is  usually  measured  by  a standard  wake 
rake  before  the  actual  design  of  the  full-scale  propeller  is  undertaken.  The  radial  distribution 
of  effective  inflow  velocity,  an  important  input  in  the  design  of  a wake-adapted  propeller,  is 
often  scaled  up  or  down  from  the  measured  circumferential-mean  velocity  profile  by  a 
constant  factor.  This  factor  is  sometimes  taken  to  be  the  ratio  of  the  Taylor  wake  fraction 
(determined  from  the  results  of  behind-condition  propulsion  experiments  on  the  basis  of 
thrust  identity  with  propeller  open-water  data  to  the  measured  volume-mean  nominal  wake 
fraction).  This  constant-factor  empirical  approach  can  be  used  only  as  a rough  approximation. 
Since  few  velocity  surveys  have  ever  been  made  in  the  presence  of  a propeller,  the  true 
distribution  of  the  effective  inflow  velocity  to  the  propeller  has  been  largely  unknown.  The 
effective  inflow  velocity  will  no  doubt  affect  propeller  powering  performance,  cavitation 
characteristics,  and  unsteady  forces  and  moments.  It  is  thus  important  that  flow  tields  in  the 
presence  of  a propeller  be  investigated  over  a wide  range  of  stern  boundary-layer 
characteristics. 

Thrust  deduction  has  received  the  attention  of  many  investigators,1 2 3 4'5  and  the  theoretical, 
experimental,  and  empirical  literature  on  the  subject  is  very  extensive;  the  many  references 
cited  by  Nowacki  and  Sharma6  have  application  to  surface  ships. 


1 Beveridge,  J.L.,  “Analytical  Prediction  of  Thrust  Deduction  for  Submersibles  and  Surface  Ships,”  Journal 
of  Ship  Research,  Vol.  13,  No.  4,  pp.  258-271  (1969). 

2Weinblum.  G.P..  “The  Thrust  Deduction,”  Journal  of  American  Society  of  Naval  Engineers,  Vol.  63, 
pp.  363-380  (1951). 

3Amtsberg,  H.,  “Investigations  on  the  Interaction  between  Hull  and  Propeller  of  Bodies  of  Revolution,” 

(in  German),  Jahrbuch  der  Schiffbautechnischen  Gesselschaft.  Vol.  54,  pp.  117-152  (I960);  also  available 
as  David  Taylor  Model  Basin  Translation  309  (1965). 

4Tsak  onas,  S.  and  W.R.  Jacobs,  “Potential  and  Viscous  Parts  of  the  Thrust  Deduction  and  Wake  Fraction 
for  an  Ellipsoid  of  Revolution,”  Journal  of  Ship  Research.  Vol.  4,  No.  2,  pp.  1-16  (I960). 

5 Nowacki,  H.,  “Potential  Wake  and  Thrust  Deduction  Calculations  for  Ship-Like  Bodies.”  Jahrbuch  der 
Schiffbautechnischen  Gesselschaft.  Vol.  57.  pp.  330-363  (1963). 

6 Nowacki.  H.  and  S.D.  Sharma.  “Surface  Effects  in  Hull  Propeller  Interaction.”  The  Ninth  Office  of  Naval 
Research  Symposium  on  Naval  Hydrodynamics,  Paris.  France  ( Aug  1972);  available  in  U.S.  Government 
Printing  Office.  ACR-203,  Vol.  2.  pp.  1845-1961  (1972). 
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In  contrast,  literature  on  the  interaction  between  a propeller  and  stern  boundary  layer  is 
limited  to  isolated  experimental  and  empirical  data.7,8 

The  first  known  theoretical  effort  to  compute  the  effective  wake  distribution  from 
measured  values  of  the  nominal  wake  and  static  pressure  distributions  was  made  in  1972  by 
D.M.  Nelson  of  the  Naval  Undersea  Center,  San  Diego.*  In  his  unpublished  work,  Nelson 
integrated  the  effective  wake  calculation  into  a propeller-design  computer  program  which  has 
been  extensively  used  to  design  propellers  for  small  submerged  bodies.  Although  the  approach 
used  in  the  present  report  to  treat  the  theoretical  derivation  of  propeller/hull  interaction  is 
independent  of  that  used  by  Nelson,  his  numerical  results  did  help  to  motivate  further  work 
on  the  problem.  Of  greatest  importance,  the  present  work  represents  the  first  time  that 
theoretical  wake  predictions  have  been  correlated  with  measured  wake  distributions  in  the 
presence  of  a propeller. 

Axisymmetrical  bodies  were  chosen  for  this  investigation  in  order  to  focus  on  the 
physical  nature  of  the  complex  interaction  between  a propeller  and  a thick  stern  boundary 
layer.  Their  geometric  simplicity  offers  considerable  experimental  and  computational 
convenience  in  treating  the  fundamental  aspects  of  the  interaction.  Definitive  experimental 
results  on  boundary-layer  characteristics  and  pressure  distributions  for  axisymmetrical  after- 
bodies with  and  without  propellers  in  operation  are  crucial  to  a better  understanding  of  the 
propeller-stern  interaction.  Thus,  much  effort  and  care  were  given  to  the  experimental 
equipment  and  procedures  described  here.  The  Laser  Doppler  Velocimeter  was  successfully 
used  to  measure  velocity  profiles  very  close  to  the  propeller.  The  measured  differences 
between  the  velocity  profiles  with  and  without  an  operating  propeller  provide  the  necessary 
clues  to  a proper  understanding  of  the  propeller/stern  boundary-layer  interaction. 

The  experimental  techniques  are  described  in  detail  and  the  results  of  the  experiments 
are  then  used  to  evaluate  several  new  analytical  techniques  as  well  as  an  existing  one.  The 
first  of  the  new  analytical  techniques  enables  the  pressure  distribution,  shear  stress  distribution, 
and  velocity  profiles  to  be  determined  over  the  entire  body  by  considering  the  interaction 
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between  the  boundary-layer  flow  next  to  the  body  and  the  external  potential  flow  in  the 
absence  of  a propeller.  The  various  components  of  the  technique  are  described  in  detail 


7Hucho,  W.-H.,  “On  the  Effect  of  a Stern  Propeller  on  the  Pressure  Distribution  and  the  Boundary-Layer  of 
an  Afterbody  of  Revolution,”  Institut  fur  Stromungsmechanik  der  Technischen  Hochschule  Braunschweign. 
Bericht  64/45,  Pt  II  (1965). 
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together  with  an  iteration  scheme  for  computing  the  interaction.  The  second  new  technique 
solves  for  the  influence  of  the  propeller  on  the  nominal  velocity  distribution  at  the  propeller 
disk  plane  in  the  absence  of  the  propeller.  The  various  assumptions  used  are  clearly  stated. 

One  of  the  principal  assumptions  is  that  the  interaction  between  the  propeller  and  stem 
boundary'  layer  may  be  taken  to  be  inviscid.  An  existing  analytical  technique  is  evaluated  for 
the  calculation  of  thrust  deduction  where  potential  flow  methods  are  used  to  calculate  the 
propeller-hull  interaction. 

MODELS  AND  EXPERIMENTAL  METHODS 
BODY  OF  REVOLUTION  MODELS 

After  preliminary  computer  evaluation  of  the  flow  characteristics  of  bodies  of  revolution, 
three  axisymmetric  afterbodies  were  selected  for  the  present  experimental  investigation.  Their 
afterbody  length/diameter  ratios  (LA/D)  were  4.308,  2.247,  and  1.484;  corresponding 
prismatic  ratios  (CpA)  were  respectively  0.606,  0.526,  and  0.416.  As  shown  in  Figure  1 and 
Table  1,  each  afterbody  was  connected  to  a parallel  middle  body  of  length  LM  and  an  existing 
streamlined  forebody  with  a bow-entrance  length/diameter  ratio  (l^/D)  of  1.82.  The  total 
length  of  the  model  was  fixed  at  a constant  value  of  3.066  m.  Other  hull  particulars  of  the 
three  models  are  listed  in  Table  1. 

The  common  forebody  and  a portion  of  the  parallel  middle  body  were  constructed  of 
wood.  The  afterbody  and  the  remaining  portions  of  the  parallel  middle  body  were  constructed 
of  molded  fiberglass;  specified  profile  tolerances  were  held  to  less  than  ±0.4  mm,  all 
imperfections  were  removed,  meridians  were  faired,  and  the  fiberglass  was  polished  to  a 
0.64-micron  rms  surface  finish. 

As  can  be  seen  from  Figure  1,  the  fullnesses  of  the  three  sterns  were  selected  to  provide 
large  variations  of  stern  flow.  Afterbody  3 (LA/D  = 1.4841  was  the  bluntest  and  had  a cosine 
curve  for  the  dimensionless  body  radius  rQ/ L of  a form  similar  to  that  used  by  Kempf.9 

ro  [*0.977137  - X/L  “I 

— = 0.025886  - 0.0196736  cos  ^ — - - p-- — ttJ 

for  0.86481 1 < X/L  < 0.977137.  This  afterbody  had  an  inflection  point  at  X/L  = 0.920974. 


9 Kempf,  G.,  “Wirbelablosung  vei  Volligen  Schiffsformen,”  Schiff  und  Hafen,  Vol.  6.  No.  7,  pp.  407-408 
(1954). 
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Figure  1 - The  Axisymmetric  Afterbodies 


TABLE  1 - HULL  PARTICULARS 
FOR  AXISYMMETRIC  AFTERBODIES 


Model 

Afterbody  1 

le/d 

1.8182 

Lm/D 

4.8485 

la/d 

4.3078 

CPE 

0.646 

^PA 

0.606 

Cp 

0.787 

ttD2/4L2 

0.006521 

ttD2/4S 

0.02675 

Y2/3/s 

0.122 

L/D 

10.9745 

S(M2) 

2.291 

UM) 

3.066 

Afterbody 


1.8182 

6.9091 

2.2472 

0.646 

0.526 

0.844 

0.006521 

0.02560 

0.122 

10.9745 

2.395 

3.066 


Afterbody  3 


0.006521 

0.02540 


10.9745 


[ 


The  offsets  of  the  two  remaining  afterbodies  are  described  in  terms  of  the  Granville 
family  of  polynomials10  by  the  expression 

t?2  =S,2[£2(£-  l)4]  +K,[|  €3  (£  - l)3]  + 1 -((-  1 )4  (10  $2  +4£  + 1) 

where  tj  = 2rQ/D  and  £ = (L  - X)/LA-  The  selected  values  of  the  coefficients  were  S,2  = 7.9254 
and  3.4,  K,  = 17.281  and  32.5,  and  LA'D/  = 4.242  and  1.812  for  Afterbodies  1 and  2, 
respectively.  The  tail  ends  of  the  three  afterbodies  were  modified  slightly  to  accommodate 
the  hub  of  an  existing  propeller  (Figure  1).  The  values  for  LA/D  given  in  Table  1 are  the 
final  values  corresponding  to  the  modified  afterbodies. 

The  Hess-Smith11  potential-flow  computer  program  for  an  unbounded  incompressible 
fluid  was  used  to  calculate  the  distributions  of  pressure  coefficient  C =(p-  p0)/(l/2)pU02 
on  the  three  afterbodies;  see  Figure  2.  As  would  be  expected,  the  magnitudes  of  the  stern 
negative  pressure  troughs  and  the  subsequent  adverse  pressure  gradients  increased  dramatically 
with  increasing  afterbody  fullness,  i.e.,  decreasing  LA/D.  Preliminary  boundary-layer  calcula- 
tions12,13 used  as  input  the  potential  flow  velocity  up  to  X/L  = 0.95  and  linear  extrapolation 
of  potential  flow  velocity  between  0.95  < X/L  < 1;  these  indicated  that  the  three  afterbodies 
would  provide  three  different  types  of  stern  flow.  It  was  predicted  that  Afterbody  3 would 
have  strong  shoulder  separation  at  about  X/L  = 0.92.  Afterbody  2 would  have  incipient 
separation  or  nearly  separated  flow  further  aft.  and  Afterbody  1 would  not  have  separation. 

In  a more  refined  calculation  of  the  thick  stern  boundary-layer  flow,  the  potential  flow  velocities, 
(based  on  pressure  distributions  iterated  from  consideration  of  the  boundary-layer  and  wake 
displacement  thicknesses)  indicated  that  Afterbody  2 would  not  have  separation.*  Thus,  the 
three  afterbodies  provided  a wide  range  of  stern  flow  for  the  investigation  of  the  propeller/ 
stern  boundary-layer  interaction. 


10Granville,  P.S.,  "Geometrical  Characteristics  of  Noses  and  Tails  for  Parallel  Middle  Bodies,”  NSRDC 
Report  3763  (1972). 

11  Hess,  J.L.  and  A.M.O.  Smith.  “Calculation  of  Potential  Flow  about  Arbitrary  Bodies.”  in  “Progress  in 
Aeronautical  Sciences,”  Vol.  8,  Pergamon  Press,  New  York  (1966),  Chapter  I. 

12Kerney.  K.P.  and  N.M.  White.  “Description  and  Evaluation  of  a Digital-Computer  Program  for  Calculating 
the  Viscous  Drag  of  Bodies  of  Revolution,”  DTNSRDC  Report  4641  (1975). 

13Cebeci,  T.  and  A.M.O.  Smith.  “Analysis  of  Turbulent  Boundary  Layers,”  Academic  Press,  New  York 
(1974),  pp.  329-384. 

•Experimental  results,  which  will  be  discussed  in  detail  later,  confirmed  that  Afterbody  3 had  shoulder 

separation  at  X/L  *0.918  and  that  Afterbodies  1 and  2 had  no  separation. 


Figure  2 - Computed  Distributions  of  Pressure  Coefficient 
on  the  Axisymmetric  Afterbodies 


MODEL  PROPELLER 

An  existing  seven-bladcd  propeller,  with  a diameter  of  54.5  percent  of  the  hull  diameter, 
was  located  at  X/L  = 0.983  for  the  experimental  investigation  of  propeller  and  stern  boundary- 
layer  interaction.  Table  2 summarizes  the  propeller  geometry.  This  propeller  was  designed 
for  a wake  distribution  which  was  different  from  the  wakes  generated  by  the  three  afterbodies. 
Thus,  an  inverse  propeller  computing  program  was  required  to  calculate  the  hydrodynamic 
characteristics  of  the  propeller  with  the  given  propeller  geometry  and  the  specified  wake 
distribution.  This  calculation  procedure  will  be  given  later  in  the  report. 
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TABLE  2 - PROPELLER  GEOMETRY 

Number  of  Blades  7 

Expanded  Area  Ratio  0.584 

Section  Meanline  NACA  a = 0.8 

Section  Thickness  Distribution  NACA  66  (with  DTNSRDC  modified  nose  and  tail) 

Rake  Angle,  deg  6.964 

Skew,  deg  30 


r/Rp 

C/Dp 

t/c 

P/Dp 

fM/c 

0.2106 

0.171 

0.235 

0.823 

0.0014 

0.3 

0.177 

0.209 

0.980 

0.0175 

0.4 

0.182 

0.182 

1.151 

0.0288 

0.5 

0.185 

0.158 

1.243 

0.0337 

0.6 

0.185 

0.135 

1.264 

0.0341 

0.7 

0.180 

0.1 16 

1.248 

0.031 1 

0.8 

0.164 

0.0995 

1.206 

0.0246 

0.9 

0.132 

0.0875 

1.157 

0.0141 

1.00 

0.069 

0.0813 

1.108 

0. 

EXPERIMENTAL  TECHNIQUE 

The  experimental  investigation  was  conducted  in  the  wind  tunnel  of  the  DTNSRDC 
anechoic  flow  facility.  The  wind  tunnel  has  a 2.438-  by  2.438-m  test  section  with  a maximum 
air  speed  of  61  m/s.  Air  is  circulated  around  a closed  loop  and  passes  through  a closed-jet 
test  section  into  an  open-jet  anechoic  chamber.  With  a model  in  place,  the  free-stream 
turbulence  level  in  the  test  section  was  measured  by  a D1SA  55M  series,  constant-temperature 
anemometer  with  linearizer.  The  diameter  of  the  tungsten  wire  used  was  0.0006  mm.  The 
measured  turbulence  levels  (^/ U'2  /U0 ) x 100  were  0.075,  0.090.  0.100.  and  from  0.12  to 
0.15  for  free-stream  velocities  UQ  of  24.4,  30.5,  38.1,  and  45.7  m/s.  respectively.  Integration 
of  the  measured  noise  spectrum  levels  in  the  test  section  from  0 to  10,000  Hz  indicated  that 
the  typical  background  acoustic  noise  at  30.5  m/s  was  around  93  dB  re  0.0002  dyn/cm2. 

These  levels  of  ambient  turbulence  and  acoustic  noises  were  considered  low  enough  so  as  not 
to  unfavorably  affect  the  measurements  of  boundary-layer  characteristics. 

The  model  was  supported  by  two  streamlined  struts  separated  by  roughly  one-third  of 
the  model  length.  The  upstream  strut  had  a 1 5-cm  chord  and  the  downstream  strut  a 3-cm 
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chord.  The  disturbances  generated  by  the  supporting  struts  were  within  the  region  below  the 
horizontal  centerplane.  Prior  to  the  experiment,  pressure  taps  and  Preston  tubes  were  used  to 
check  the  axisymmetrical  characteristics  of  the  stern  flow  at  X/L  = 0.90,  0.95,  and  0.977. 

The  circumferential  variations  of  pressure  and  surface  shear  stress  on  the  upper  half  of  the 
body  at  these  three  locations  were  within  two  percent.  All  the  measurements  were  made  in 
the  vertical  centerplane  along  the  upper  meridian  and  there  was  little  extraneous  effect  from 
the  supporting  strut.  Each  stern  prou  3.  from  the  closed-jet  working  section  of  the  wind 
tunnel  into  the  anechoic  chamber  (6.4  x 6.4  x 6.4  m)  located  upstream  of  the  diffuser.  The 
propeller  was  driven  by  a 9-kW  high-speed  motor  mounted  inside  the  stern  of  the  model.  The 
rate  of  revolution  of  the  propeller  shaft  was  measured  by  a magnetic  pickup  and  displayed  by 
a digital  counter. 

Measurements  of  surface  shear  stress,  surface  pressure  distribution,  and  boundary-layer 
axial  velocity  distribution  were  made  on  each  stern  first  without  a propeller  and  then  with  a 
propeller  operating  at  one  or  two  advance  coefficients.  The  following  equipment  was  used. 

A 1.83-mm  Preston  tube  was  taped  to  the  sterns  at  successively  further  aft  locations  in  order 
to  measure  shear  stress  at  1 2 locations  along  the  upper  meridian  of  each  stern.  (The  tube  has 
been  calibrated  in  a 2.54-cm  water  pipe  flow  facility  described  by  Huang  and  von  Kerczek.14) 
Because  a logarithmic  boundary-layer  velocity  profile  is  required  in  order  to  analyze  Preston 
tube  data,  shear  stresses  could  not  be  accurately  determined  by  this  method  at  locations 
downstream  of  separation.  Pressure  taps  (0.8-mm  diameter)  were  used  to  measure  steady 
pressure  at  the  same  locations  as  the  Preston  tubes.  The  holes  were  connected  by  “Tygon” 
plastic  tubing  to  a scanning  valve  and  then  to  a strain-gage  type  pressure  transducer  located 
underneath  the  floor  of  the  closed-jet  test  section.  The  pressure  transducer  (ACEC  Model  310) 
was  used  in  conjunction  with  a signal  conditioner  (Endevco  4470).  an  amplifier  (Dana  Model 
2856),  and  a digital  voltmeter  to  measure  the  pressure  from  Preston  tubes  and  surface  pressure 
taps.  An  automatic  control  for  the  scanning  valve  was  located  outside  the  wind  tunnel. 

Steady  boundary-layer  axial  velocity  profiles  were  measured  by  a laser  Doppler  velocimeter 
(LDV)  at  five  stations  on  each  stern.  The  LDV  was  located  on  an  optical  bench  in  the 
quiescent  region  of  the  anechoic  chamber  and  was  operated  in  a dual-beam  off-axis  backscatter 
mode,  which  is  described  in  detail  in  the  appendix.  Figure  3 shows  photographic  views  of  the 
optical  arrangements  for  two  of  the  models. 


l4IIuang.  T.T.  and  C.H.  von  Kcrczek.  “Shear  Stress  and  Pressure  Distribution  on  a Surface  Ship  Model: 

Theory  and  Experiment.”  The  Ninth  Office  of  Naval  Research  Symposium  on  Naval  Hydrodynamics,  Paris. 
France  (Aug  1972);  available  in  U.S.  Government  Printing  Office,  ACR-203,  Vol.  2,  pp.  1963  - 2010  ( 1972). 


The  focal  length  of  the  LDV  optics  was  1.5  m and  the  total  beam  angle  was  3.722  deg.  The 
effective  probe  volume  of  this  off-axis  backscatter  optics  was  roughly  ellipsoidal,  with 
dimensions  of  about  0.5  x 5 mm,  and  the  probe  volume  could  be  focused  in  the  stern  boundary 
layer  at  distances  down  to  2 mm  from  the  hull  surface. 

A detailed  analysis  of  the  accuracy  of  the  measurements  has  not  been  made.  However, 
the  standard  deviations  of  the  measured  data  can  be  estimated  from  repeat  runs.  The  standard 
deviations  of  the  measured  static  wall  pressure  and  shear  stress  were  about  5 percent  of  their 
mean  values  and  the  standard  deviations  of  the  measured  velocities  were  about  2 percent  of 
the  free-stream  velocity. 


METHOD  FOR  CALCULATING 
POTENTIAL-FLOW-BOUNDARY-LAYER  INTERACTION 

SURVEY  OF  PREVIOUS  METHODS 

Several  studies  have  recently  attempted  to  calculate  the  complete  pressure  field  acting  on 
axisymmetric  bodies,  including  the  effects  of  viscosity.  They  generally  start  with  an  initial 
estimate  of  the  pressure  distribution,  usually  the  potential  flow  distribution  with  some 
modifications  in  the  tail  region,  and  employ  standard  boundary-layer  theory  to  calculate  the 
flow  over  the  forward  region  of  the  body.  The  boundary-layer  equations  may  be  in  either 
integral  or  differential  form.  These  methods  differ  greatly,  however,  in  the  approach  used  to 
calculate  the  flow  field  at  the  stern  of  the  body  (where  the  boundary  layer  may  be  thick),  in 
the  near  wake  region,  and  in  the  far  wake  region.  These  differences  will  be  pointed  out  below. 
The  initial  flow  calculations  are  then  used  to  modify  the  geometry  of  the  body  and  wake, 
usually  by  adding  the  local  displacement  thickness  as  suggested  by  Lighthill. 15  Potential-flow 
methods  are  then  utilized  to  compute  the  pressure  distribution  around  the  modified  body. 

Flow  calculations  are  repeated  for  the  new  pressure  distribution,  and  this  basic  scheme  is 
continued  until  the  pressure  distributions  from  two  successive  approximations  agree  to  within 
a given  error  criterion. 

The  modeling  of  the  flow  over  the  stern  of  the  body  and  in  the  wake  region  ranges  from 
very  simple  to  very  complex.  Perhaps  the  simplest  modeling  is  by  Beatty.16  He  merely  adds 
the  displacement  thickness  calculated  by  the  Douglas  CS  differential  boundary-layer 


15  Lighthill,  M.J.,  “On  Displacement  Thickness,”  Journal  of  Fluid  Mechanics.  Vol.  4.  Pt  4.  pp.  383-392  (1958). 

16 Beatty,  T.D.,  “A  Theoretical  Method  for  the  Analysis  and  Design  of  Axisymmetric  Bodies.”  National 
Aeronautics  and  Space  Administration  CR-2498  (1975). 
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method13,17  to  the  body  up  to  the  trailing  edge.  He  then  calculates  the  potential  flow 
distribution  about  the  resulting  open-ended  body  and  a semi-infinite  stream  tube  (which  leaves 
the  trailing  edge  of  the  open  body  and  proceeds  downstream)  until  the  cross  section  of  the 
stream  tube  approaches  a constant  value  far  downstream.  In  the  case  of  separated  flow,  the 
portion  of  the  displaced  body  from  the  separation  point  to  the  trailing  edge  is  assumed  to  be 
an  open-ended  cylinder  with  radius  equal  to  the  radius  of  the  original  body  plus  the  displace- 
ment thickness  at  the  separation  point.  A circular  arc  is  used  to  fair  the  displaced  body 
ahead  of  the  separation  point  into  the  cylindrical  body.  Agreement  is  good  with  experimental 
results  for  the  case  of  a streamlined  body  with  fineness  ratio  12:1,  but  it  is  only  fair  for  the 
case  of  a sphere  at  a subcritical  Reynolds  number. 

Beveridge18  considers  a generally  similar  model  but  allows  the  wake  to  extend  to  infinity. 
However,  he  does  not  specify  such  details  as  how  the  body  is  faired  into  the  wake  and  the 
width  of  the  wake  at  infinity.  Agreement  of  this  method  is  relatively  poor  with  experimental 
results  for  a Series  58  body  with  fineness  ratio  10. 

In  a general  method  outlined  by  Cebeci,  Mosinskis,  and  Smith17  for  adding  the  displace- 
ment thickness  to  the  original  body,  the  modeling  of  the  wake  is  left  unspecified.  The 
difference  between  the  potential  flow  and  viscous  pressure  distributions  is  shown  for  a few 
cases,  but  no  comparisons  with  experiment  are  given. 

The  Myring19  method  is  somewhat  more  complex.  He  uses  boundary-layer  equations  in 
integral  form  to  solve  for  the  flow  on  the  body  and  in  the  wake.  The  principal  difference 
between  the  equations  for  the  two  flow  regions  lies  in  the  form  for  the  entrainment  coefficient 
which  is  larger  in  the  far  wake  because  of  increased  turbulent  mixing.  In  the  region  close  to 
the  tail  of  the  body,  the  entrainment  coefficient  appropriate  to  the  body  is  faired  into  the 
corresponding  coefficient  for  the  far  wake  by  means  of  an  exponential  function.  However,  an 
important  constant  in  the  fairing  function  is  unspecified  and  must  be  determined  by  comparison 
with  available  experimental  data. 


17Cebeci,  T.  et  al.,  “Calculation  of  Viscous  Diag  and  Turbulent  Boundary-Layer  Separation  on  Two- 
Dimensional  and  Axisymmetric  Bodies  in  Incompressible  Flows.”  Douglas  Aircraft  Report  MDC  J0973-01 
(1970). 

18Beveridge,  J.L.,  “Pressure  Distribution  on  Towed  and  Propelled  Streamlined  Bodies  of  Revolution  at  Deep 
Submergence,”  David  Taylor  Model  Basin  Report  1665  (1966). 

19Myring,  D.F.,  “The  Profile  Drag  of  Bodies  of  Revolution  in  Subsonic  Axisymmetric  Flow.”  Royal  Aircraft 
Establishment  (Great  Britain)  Technical  Report  72234  (1972). 
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The  computed  pressure  distributions  agree  well  with  experimental  results  foi  the  Lyon 
body20,21  and  the  airship  AKRON.22 

Nakayama23  considers  an  even  more  complex  model.  He  uses  standard  boundary-layer 
equations  in  integral  form  to  model  the  How  over  the  forward  portion  of  the  body  and  in 
the  far  wake.  A significant  difference  between  this  model  and  those  previously  mentioned  is 
that  in  the  region  near  the  tail  and  in  the  near  wake,  allowance  is  made  for  the  transverse 
variation  of  pressure  across  the  boundary  layer.  Also,  a simple  linear  profile  is  assumed  for 
the  normal  velocity.  Another  significant  difference  is  that  the  displacement  thickness,  which 
is  an  integrated  effect  of  the  boundary  layer,  is  not  used  to  modify  the  body.  Instead,  the 
potential  solution  is  matched  to  the  boundary  layer  calculations  at  the  edge  of  the  boundary 
layer  and  wake.  The  computed  pressure  distributions  are  in  good  agreement  with  experimental 
results  for  the  case  of  the  Lyon  body,20-21  but  in  only  fair  agreement  for  the  case  of  a 
modified  spheroid.24 

Perhaps  the  most  complex  method  outside  of  a complete  solution  of  the  Navier-Stokes 
equations  in  the  entire  flow  field  is  that  considered  by  Brune,  Rubbert,  and  Forester.25  This 
method  may  be  best  described  by  referring  to  their  sample  case  of  a 4:1  prolate  ellipsoid. 
Standard  boundary-layer  equations  in  differential  form  are  used  to  calculate  the  How  over  the 
front  half  of  the  body,  0 < X/L  < 0.5,  where  X is  the  axial  distance  measured  from  the  ni  *e 
and  L is  the  total  length  of  the  body.  The  complete  Navier-Stokes  equations  were  used  to 
calculate  the  flow  over  the  rear  half  of  the  body  and  in  the  near  wake.  0.5  < X/L  < 1.26. 

The  shape  of  the  displacement  tail  for  X/L>  1.26  is  taken  to  be  a faired  shape  which  ends  in 
a point  at  X/L  = 2.5.  The  authors  give  no  precise  method  for  obtaining  this  displacement  tail 
but  simply  state  that  it  is  specified  by  “extrapolation.”  No  comparisons  with  experimental 
results  are  given.  It  is  worthwhile  to  note  that  use  of  the  Navier-Stokes  equations  greatly 


20Lyon,  H.M..  “Effect  of  Turbulence  on  Drag  of  Airship  Models,”  Aeronautical  Research  Committee 
(Great  Britain)  Reports  and  Memoranda  1511  (1932). 

7 I 

Lyon.  H.M.,  “A  Study  of  the  Flow  in  the  Boundary  Layer  of  Streamlined  Bodies.”  Aeronautical  Research 
Committee  Reports  and  Memoranda  1622  (1934). 

22Freeman,  H.B.,  “Pressure  Distribution  Measurements  on  the  Hull  and  Fins  of  a 1 /40-Scale  Model  of  the 
US  Airship  ‘Akron’,”  National  Advisory  Committee  for  Aeronautics  Report  443  (1932). 

23Nakayama,  A.,  “Viscid-Inviscid  Interaction  Due  to  the  Thick  Boundary  Layer  near  the  Tail  of  a Body  of 
Revolution,”  Ph.D.  Thesis,  University  of  Iowa  (1974). 

24Patel,  V.C.  et  a).,  “Measurements  in  the  Thick  Axisymmetric  Turbulent  Boundary  Liver  near  the  Tail  of  a 
Body  of  Revolution,”  Journal  of  Fluid  Mechanics,  Vol.  63,  Pt  2.  pp.  345-367  (1974). 

25  Brune,  G.W.  et  al.,  “The  Analysis  of  Flow  Fields  with  Separation  by  Numerical  Matching.”  Symposium  on 
Flow  Separation,  Advisory  Group  for  Aerospace  Research  and  Development.  Germany  (1975). 
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increases  the  overall  computer  time  requirements  of  the  method.  The  authors  state  that  the 
laminar  boundary  layer  and  potential  flow  calculations  require  about  40  sec  of  CDC-6600  time 
per  iteration  whereas  the  Navier-Stokes  equations  require  approximately  40  min  per  iteration. 
This  method  is  limited  to  laminar  flow.  The  extension  to  turbulent  flow  would,  of  course, 
require  even  more  computer  time. 


PRESENT  ITERATION  SCHEME 
Outline  of  the  Present  Scheme 

The  Douglas  CS  differential13, 17  boundary-layer  method  is  used  to  calculate  the  How 
over  the  body,  and  the  integral  wake  relations  given  by  Granville26  are  used  to  calculate  the 
flow  in  the  wake.  The  calculated  displacement  thicknesses  are  then  utilized  to  generate  an 
overall  body-wake  displacement  model.  However,  since  neither  of  the  two  methods  properly 
models  the  thick  boundary  layer  in  the  stern/near-wake  region,  where  there  may  be  sizeable 
pressure  variations  across  the  boundary  layer,  the  calculated  displacement  body  is  not  assumed 
to  be  valid  in  the  region  0.95  < X/L  < 1.05.  Instead,  a fifth-degree  polynomial  is  used  in 
this  region,  with  the  six  constants  determined  by  requiring  the  thickness,  slope,  and  curvature 
to  be  equal  to  those  calculated  at  X/L  = 0.95  and  1.05.  Should  separation  occur  before 
X/L  = 0.95,  the  upstream  matching  point  is  moved  to  the  separation  point.  The  pressure 
distribution  over  the  displacement  body  is  calculated  by  the  Douglas-Neumann  method.* 1,27 
The  resulting  pressure  distribution  is  then  used  to  calculate  the  new  viscous  flow  over  the 
body  and  in  the  wake.  This  process  is  repeated  until  two  successive  pressure  distributions 
agree  to  within  a specified  error  criterion  everywhere  in  the  flow  field. 

The  potential-flow  pressure  distribution  over  the  body,  without  a displacement  thickness 
correction,  constitutes  the  first  iteration  except  that  the  velocity  ratio  is  linearly  extrapolated 
over  the  last  5 percent  of  body  length.  In  the  wake,  the  velocity  ratio  is  initially  taken  to 
increase  exponentially  from  the  value  at  the  tail  to  1 in  the  far  wake. 

Since  the  Douglas-Neumann  method  and  the  Douglas  CS  differential  boundary-layer  method 
are  widely  used  to  calculate  flows  around  axisymmetric  bodies,  only  brief  descriptions  are  given 
here.  Other  aspects  of  the  method  are  described  in  greater  detail. 

26Granville.  P.S.,  “The  Calculation  of  the  Viscous  Drag  of  Bodies  of  Revolution,”  David  Taylor  Model  Basin 
Report  849  (1953). 

27Smith,  A.M.O.  and  J.  Pierce.  “Exact  Solution  of  the  Neumann  Problem.  Calculation  of  Non-Circulatory 
Plane  and  Axially  Symmetrical  Flows  About  or  Within  Arbitrary  Boundaries.”  Douglas  Aircraft  Report 
F.S-26988  (1958). 
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Douglas-Neumann  Computer  Program 


In  this  method, 11,27  the  Laplace  equation  for  the  unknown  potential  is  converted  into  a 
Fredholm  integral  equation  of  the  second  kind  for  the  unknown  source  density  distribution 
on  the  surface  of  the  body.  The  body  is  considered  as  composed  of  a series  ot  frustums  ot 
cones,  and  the  integral  equation  is  converted  to  a set  of  algebraic  equations  for  the  unknown 
constant  source  density  on  the  surface  of  each  frustrum.  Once  the  source  densities  have  been 
obtained,  it  is  a simple  matter  to  obtain  all  flow  variables  of  interest. 

The  original  computer  program  included  several  options  for  solving  the  algebraic  equations. 
For  the  number  of  points  typically  used  in  the  present  calculations  (about  140  on  the  body 
and  35  in  the  wake),  the  Seidel  iterative  procedure  is  the  most  economical. 


Douglas  CS  Differential  Boundary- Layer  Method 


This  method12,17  solves  the  boundary-layer  equations  in  partial  differential  equation  form. 
Thus,  no  assumptions  are  made  about  the  form  of  the  velocity  profiles.  The  method  also 
accounts  for  the  effects  of  transverse  curvature,  which  are  important  in  the  tail  region  where 
the  thickness  of  the  boundary  layer  may  be  greater  than  the  body  radius.  However,  the 
pressure  variation  across  the  boundary  layer,  which  may  be  appreciable  in  the  tail  region,  is 
neglected.  The  equations  are  solved  by  the  Keller  box  method12  which  essentially  consists  ot 
first  reducing  the  original  set  of  second-order  partial  differential  equations  to  a tirst-order 
system  and  then  solving  the  resulting  set  by  an  implicit  finite  difference  method. 

The  location  of  transition  from  laminar  to  turbulent  flow  may  be  specified.  In  the 
experimental  work  described  here,  a 0.61-mm-diameter  trip  wire  was  used  at  an  X/L  location 
of  5 percent.  The  Granville26  boundary-layer  computation  method  programmed  by  Kerney 
and  White12  is  used  to  compute  the  laminar  boundary  layer  upstream  ol  the  trip  wire  and  the 
turbulent  boundary  layer  downstream  of  the  trip  wire.  The  jump  in  momentum  thickness 
Reynolds  number  ARy  due  to  the  parasitic  drag  of  the  wire  is  specified  according  to 
McCarthy,  Power,  and  Huang28  by 


ukK 


^McCarthy.  J.H.  et  al„  “The  Roles  of  Transition,  Laminar  Separation,  and  Turbulence  Stimulation  in  the 
Analysis  of  Axisymmetric  Body  Dtag."  The  I levcnth  Office  of  Naval  Research  Symposium  on  Naval 
Hydrodynamics,  London,  Fngland  (ll)76). 


•• 

I 


15 


Here  uK  is  the  local  velocity  at  the  tip  of  the  wire.  Ue  is  the  potential  flow  velocity  at  the 
trip  location,  K is  the  diameter  of  the  wire  (0.61  mm),  and  Cp  is  the  drag  coefficient  of  the 
wire  (taken  as  0.75). 

When  the  flow  was  probed  with  a hot  ware,  the  trip  wire  was  found  to  effectively  trip 
the  flow  at  1 cm  downstream  from  the  wire.  Because  of  the  finite  parasitic  drag  of  the  wire, 
the  boundary  layer  can  be  considered  to  become  turbulent  at  a virtual  origin  upstream  of  the 
trip  wire.  This  virtual  origin  for  the  turbulent  flow  is  defined  such  that  the  sum  of  the 
laminar  frictional  drag  from  the  nose  to  the  trip,  the  parasitic  drag  ot  the  trip,  and  the  turbu- 
lent frictional  drag  after  the  trip  equals  the  sum  of  the  laminar  frictional  drag  from  the  nose 
to  the  virtual  origin  and  the  turbulent  frictional  drag  from  the  virtual  origin  to  the  after  end 
of  the  model.28  As  can  be  seen  from  Figure  4,  the  virtual  origin  can  be  obtained  by  finding 
the  intersection  of  lines  L and  T (plotted  as  r0^/rniax  versus  X/L,  where  rQ  is  the  model 
local  radius  and  r v is  the  maximum  radius  of  the  body).  Line  L is  valid  for  the  axisymme- 
trical  laminar  boundary  layer  and  line  T is  valid  for  the  axisymmetrical  turbulent  boundary 
layer  with  a jump  of  ARy  above  Line  L at  the  trip  location.  The  magnitude  of  the  jump 
ARy  is  related  to  the  parasitic  drag  of  the  trip  wire.28  The  location  of  transition  in  the 
mathematical  model  for  the  present  boundary-layer  calculation  is  specified  at  the  virtual 
origin.  The  initial  value  of  Ry  at  the  virtual  origin  for  turbulent  boundary-layer  computation 
is  then  equal  to  the  laminar  value  ol  R^  (Line  L)  at  the  virtual  origin. 


The  Granville  Wake  Relations 


The  differentia)  momentum  equation  for  wake  flow  is  simply  the  boundary-layer  equation 
with  skin  friction  neglected26 


dS7  . . S2  dLfe 

dX  + M + -’  II  dX 


0 


( 1 ) 


where 


SI  = momentum  area  = 

5 = boundary-layer  thickness 

u = tangential  component  of  the  mean  fluid  velocity  at  any  point  inside 
the  boundary  layer 

l!  = computed  potential  flow  velocity  at  the  boundary  of  the  displacement 
body  representing  the  wake 

r = radial  distance  measured  from  the  body  axis 
y = distance  normal  to  the  body  surtace 
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h = axisymmctric  shape  factor  = A */Sl 
A*  = displacement  area  = J ^1  - jj-j  r dy 


If  ft  is  considered  to  be  the  unknown,  then  h and  Ue  must  be  specified.  Ue  may  be  obtained 
from  the  calculation  of  the  pressure  distribution  for  the  displacement  body  of  the  previous 
iteration  by  using  the  Bernoulli  equation 


where  UQ  is  the  free-stream  velocity  and  Cp  is  the  pressure  coefficient.  By  analyzing  the  wake 
data  behind  a Mark  13  torpedo.29  Granville26  postulates  an  empirical  equation  relating  U£  and 
h.  This  equation  may  be  solved  to  obtain  h 


h = 1 +(h{  - 1) 


Ul  «/< 


*■«? 


where  the  subscript  t denotes  conditions  at  the  tail  of  the  body;  Granville  recommends  the 
value  q = 7. 

All  that  remains  is  to  specify  an  initial  condition  for  ft,  namely,  the  value  of  ft  at  the 
beginning  of  the  wake.  It  would  be  incorrect  to  take  this  value  as  that  computed  by  the 
Douglas  CS  method  at  the  tail  of  the  body  since,  as  mentioned  previously,  this  method  does 
not  accurately  model  the  thick  boundary  layer  at  the  stern.  Instead,  the  approach  taken  is 
first  to  compute  the  value  of  ft  in  the  far  wake  based  on  the  calculated  drag  coefficient.  This 
may  be  done  by  equating  the  total  drag  on  the  body  to  the  net  rate  of  loss  of  momentum  of 
the  flow  in  the  axial  direction.  In  view  of  the  definition  given  for  ft  in  Equation  (2),  this 
equation  may  be  expressed  in  terms  of  ftQ,  the  far  wake  momentum  area,  and  C D|  A.  the  drag 
coefficient  based  on  frontal  area26: 


2»rft0pU02  DI  A 2 ^ 


’Eggers.  H.A..  “Wake  Survey  of  the  Mark  13  Torpedo,’’  David  Taylor  Model  Basin  Report  583  (1947). 
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C r 2 

VDFA  max 


where  r ov  is  the  maximum  radius  of  the  body.  From  Equations  (1 ) and  (6),  Granville26 
max 

derives  the  following  equation  for  fit  in  terms  of  S2Q 


=/Ut\  (ht  + 2)  q + 3 


is  then  the  initial  value  of  S2  required  to  start  the  integration  of  Equation-  ( 1 ). 

Modeling  of  Displacement  Body 

The  computed  effective  displacement  thickness  a*  is  added  to  the  body  for  X/L<S  Xw/L, 
where  Xw  is  the  upstream  matching  point,  i.e.,  the  smaller  of  the  values  0.05  or  the  location 
of  separation  X$.  The  value  of  a*  is  related  to  the  displacement  area  A*  by  26 

-r  + \/ r 2 +2  A*  cos  a 

* ° v ° nm 


where  r is  the  radius  of  the  body  and  a is  the  angle  which  the  tangent  to  the  body  makes 
with  the  axis  of  the  body.  In  the  far  wake,  where  rQ  = a = 0,  a*  is  simply  given  by 


Although  the  displacement  wake  should  theoretically  be  infinite  in  length,  it  is  typically 
terminated  at  X/L  = 30  in  the  present  calculations  for  the  sake  of  computational  convenience. 
At  this  point,  the  value  of  usually  agrees  with  the  value  of  given  in  Equation  <81  to 
three  significant  figures.  The  tail  of  the  displacement  wake  is  left  open  so  that  the  tlow  leaves 
parallel  to  the  tail.  This  is  considered  preferable  to  having  a blunt  end  where  the  flow  must 
turn  around  a sharp  corner,  thus  causing  a jump  in  the  local  pressure  coefficient. 

For  X /l.  < X/L<  1.05,  where  neither  the  thin  boundary  layer  nor  far  wake  assumptions 
are  valid,  a fifth-degree  polynomial  is  used  to  represent  the  displacement  body.  The  six 
constants  of  the  polynomial  are  adjusted  to  obtain  continuity  of  displacement,  slope,  and 
curvature  at  the  upstream  and  downstream  matching  points. 


AT  X/L  =*  0.95J3 


Convergence  Behavior 

In  an  earlier  displacement  model  used  in  the  present  investigation,  only  continuity  of 
displacement  and  slope  were  enforced  at  the  matching  points.  The  computed  values  of  Cp  at 
both  matching  points  showed  large  local  jumps.  With  the  addition  of  continuity  of  curvature, 
the  pressure  distribution  for  a particular  iteration  became  smooth.  However,  it  was  found 
that  the  convergence  near  the  matching  points  was  substantially  slower  than  anywhere  else  on 
the  displacement  body.  This  is  because  the  requirements  on  slope  and  curvature,  which  are 
derivatives  of  the  displacement,  tend  to  accentuate  the  differences  between  the  displacement 
bodies  of  two  successive  iterations. 

As  typical  examples.  Figure  5 shows  the  values  of  Cp  at  the  upstream  and  downstream 
matching  points  and  Table  3 shows  the  values  of  Cp  at  other  points  on  the  body  for  each 
successive  iteration  for  Afterbody  1 at  a Reynolds  number  of  5.88  x 106. 


Figure  5 - Convergence  Characteristics  of  Boundary-Layer 
and  Potential-Flow  Iteration  Scheme 
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In  this  figure.  Iteration  n refers  to  the  pressure  distribution  for  the  nth  displacement  body 
(Iteration  0,  not  shown  in  the  figure,  would  refer  to  the  potential  flow  pressure  distribution 
over  the  original  body).  Figure  5 indicates  that  the  Cp’s  converge  slowly  to  a mean  value 
which  is  approximately  the  arithmetic  average  of  two  successive  iterations.  Even  after  eight 
iterations,  the  Cp’s  at  the  upstream  matching  point  have  not  quite  converged  to  within  a 0.01 
error.  On  the  other  hand,  Cp’s  elsewhere  in  the  flow  field  usually  converge  to  0.01  after  the 
third  iteration. 


TABLE  3 - CONVERGENCE  CHARACTERISTICS  OF  THE  BOUNDARY-LAYER 
AND  POTENTIAL  FLOW  ITERATION  SCHEME 


No.  of 

Value  of  Cp  at  Various  X/L 

Iteration 

0.9106 

0.9390 

0.9538 

0.9652 

0.974 

0.986 

1 

0.525 

0.1 123 

0.1310 

0.1234 

0.1016 

2 

0.0358 

0.0729 

0.0999 

0.1236 

0.1350 

0.1389 

3 

0.0357 

0.0977 

0.1268 

0.1366 

0.1354 

0.1225 

4 

0.0407 

0.0834 

0.1041 

0.1229 

0. 1 306 

0.1299 

5 

0.0910 

0.1280 

0.1364 

0.1375 

6 

0.0399 

0.0882 

0.1075 

0.1231 

0.1288 

7 

0.0360 

0.0874 

0.1201 

0.1361 

0.1390 

0.1302 

8 

0.0388 

0.0910 

0.1 107 

0.1233 

0. 1 275 

0.1238 

In  view  of  the  above,  most  of  the  numerical  results  presented  in  this  report  are  based  on 
three  iterations.  The  pressure  distribution  used  to  make  the  final  boundary-layer  calculations 
is  taken  to  be  the  arithmetic  average  of  the  Cp’s  from  Iterations  2 and  3. 

Calculation  of  Output  Variables 


In  addition  to  the  pressure  distribution,  several  other  variables  of  interest  are  calculated 
at  each  station  along  the  body.  These  include  most  of  the  boundary-layer  variables  calculated 
in  the  original  Douglas  CS  program.10  The  original  output  has  been  expanded  significantly  in 
two  areas. 


10Cebeci.  T.  et  al.,  “A  Finite-Difference  Method  for  Calculating  Compressible  laminar  and  Turbulent 
Boundary  Layers,”  Douglas  Aircraft  Report  DAC -67131  (1969). 
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First,  whereas  the  original  program  printed  out  only  the  tangential  velocity  profile,  the 
present  program  also  prints  out  the  normal,  axial,  and  radial  velocity  profiles.  Wake  distribu- 
tions in  the  propeller  disk  plane  are  usually  given  in  terms  of  the  latter  two  protiles. 

Second,  the  total  drag  coefficient  of  the  body  is  computed  by  three  separate  methods. 
The  first  employs  the  Squire-Young  formula  as  given  in  Cebeci  et  al.17 


4 7rfit  /Ut\(ht  + 5)/2 
Cdsy  = A 


(12) 


where  A is  the  reference  area.  The  second  method  uses  the  Granville  formula  as  given  in 
Cebeci  et  al.17 


4tt  S2t  /Ut\l7(ht  + 2)  + 31/8 

CDG  = A \lTJ 


(13) 


The  third  method  considers  the  total  drag  coefficient  as  the  sum  of  the  total  skin  friction 
coefficient  CF  and  the  pressure  drag  Cpv 

C = C + C (14) 

F pv 

where  Cf;  is  given  by17 


Here  Cj-  is  the  local  skin  friction  coefficient  = Tw/-p  Ue2  and  is  the  local  shear  stress. 

At  first,  the  pressure  drag  Cpv  was  obtained  by  integrating  the  final  viscous  pressure 
distribution  around  the  body  (usually  the  arithmetic  average  of  the  Cp’s  from  the  second  and 
third  iterations).  As  a check  on  this  method,  it  was  decided  to  also  integrate  the  original 
potential  flow  pressure  distribution  around  the  body.  It  was  tound  that  this  did  not  give  the 
theoretically  predicted  zero  pressure  drag  and  often  resulted  in  a value  of  the  same  order  of 
magnitude  as  the  viscous  pressure  drag.  A closer  inspection  of  the  detailed  output  revealed 
that  the  contributions  to  the  total  drag  from  the  first  few  stations  near  the  nose  were  of  the 
same  order  of  magnitude  as  the  final  integrated  pressure  drag.  In  view  of  this,  some  of  the 
bodies  were  rerun  with  a finer  grid  in  the  nose  region.  Also,  an  integration  method  of 
higher  accuracy.  Simpson’s  rule,  was  used  instead  of  the  simpler  trapezoidal  rule.  Although  these 
changes  usually  produced  large  changes  in  the  computed  potential  flow  drag,  the  resulting 
values  were  not  necessarily  closer  to  the  theoretically  predicted  value  ol  zero.  It  became  clear 
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that  a prohibitively  large  number  of  points  near  the  nose  would  be  required  to  obtain 
numerically  the  theoretically  predicted  value  of  zero  for  the  potential  flow  pressure  drag. 

Since  the  numerical  inaccuracy  in  the  computed  pressure  drag  was  largely  because  of  the 
pressure  distribution  near  the  nose,  which  is  essentially  not  affected  by  viscosity,  it  was 
decided  to  compute  the  viscous  pressure  drag  Cpv  by  integrating  ACpv  around  the  body,  where 


A<-pv  ^pvf  “ ^pp 


(16) 


Here  Cpvf  is  the  final  viscous  pressure  coefficient  and  Cpp  is  the  initial  potential  pressure 
coefficient.  The  pressure  drag  is  then  given  by 
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The  overall  drag  coefficient  has  been  computed  for  three  reference  areas:  frontal  area, 
wetted  area,  and  (volume)2^3.  The  frontal  area  and  (volume)2/3  are  of  interest  under  certain 
design  constraints  where  the  frontal  area  and/or  the  volume  are  taken  to  be  fixed.  On  the 
other  hand,  the  use  of  the  wetted  area  is  convenient  when  the  drag  of  a body  is  compared  to 
that  of  an  equivalent  flat  plate. 


Description  of  Computer  Program 

The  above  method  has  been  synthesized  into  a single  computer  program  by  essentially 
combining  the  Douglas  potential-flow  program27  and  the  Douglas  CS  boundary-layer 
program.13,30  Since  the  original  potential-flow  program  consisted  of  several  overlay  links,  the 
boundary-layer  program  was  made  the  last  link  of  the  potential-flow  program.  In  order  to 
simplify  the  overall  program  as  much  as  possible,  several  computations  contained  in  the 
original  potential-flow  program  were  deleted,  including  all  the  computations  for  cross  flow. 
Also,  all  of  the  options  for  solving  for  source  densities  were  deleted  except  for  the  Seidel 
iteration  method. 

The  synthesized  program  expands  the  number  of  boundary-layer  output  variables  and 
adds  the  calculation  of  flow  in  the  wake  to  the  boundary-layer  program.  The  program  has 
been  written  so  that  it  can  bypass  the  potential  flow  calculations  and  directly  accept  pressure 
distribution  as  an  input.  In  these  cases,  the  program  also  allows  for  the  reading  in  of 
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additional  data  which  modify  the  input  pressure  distribution.  This  feature  is  of  direct  interest 
for  studying  the  influence  of  the  propeller  on  the  boundary-layer  flow  over  the  body.  In 
these  cases,  the  input  pressure  distribution  may  be  the  calculated  or  measured  bare-hull 
pressure  distribution  and  the  modifying  pressure  distribution  is  that  due  to  the  propeller. 

On  the  CDC  6700  computer  currently  being  used  at  the  Center,  470  sec  of  execution 
time  are  typically  needed  to  compute  the  three  iterations  required  to  arrive  at  the  final 
pressure  distribution.  The  program  requires  approximately  65  sec  to  compile  and  a computer 
memory  of  approximately  145,000  octal  words. 

PROPELLER-HULL  INTERACTION:  FRICTIONAL  AND 
PRESSURE -DEFECT  COMPONENTS  OF  THRUST  DEDUCTION 

It  is  well  known  in  naval  architecture  that  the  delivered  propeller  thrust  is  greater  than 
the  hull  resistance  in  the  absence  of  the  propeller,  except  in  very  rare  cases.  Traditionally,  the 
increase  in  resistance  due  to  propeller-hull  interaction  has  been  defined  in  terms  of  a thrust 
deduction  fraction  t 


where  T is  the  delivered  propeller  thrust  and  R is  the  hull  resistance.  The  theoretical  and 
experimental  literature  on  thrust  deduction  is  very  extensive.1'6  Weinblum2  summarized  the 
earlier  contributions  of  hydrodynamicists  to  the  understanding  and  formulation  of  the  physics 
of  thrust  deduction.  In  addition,  an  extensive  list  of  references  with  application  to  surface 
ships  can  be  found  in  the  paper  by  Nowacki  and  Sharma.6  A potential-flow  computation 
scheme  based  on  Lagaliy’s  theorem  for  calculating  thrust  deduction  was  presented  by 
Beveridge.1 

An  operating  propeller  produces  an  upstream  suction  which  causes  a reduction  in  pressure 
on  the  afterbody  and  results  in  an  increase  of  hull  pressure  drag.  Simultaneously,  this  suction 
also  causes  an  increase  of  shear  stress  on  the  afterbody  and  results  in  an  increase  of  frictional 
drag.  Thus,  the  thrust  deduction  for  a fully  submerged,  axisymmetrical.  self-propelled  body 
may  be  divided  into  two  components,  namely,  t = tp  + tp.  The  pressure-defect  component 
tp  is  given  by 
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and  the  frictional  component  t..  by 
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where 
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(thrust  loading  coefficient  based  on  ship  speed) 


AC  = (Cn),  — (Cn ) (decrease  of  pressure  coefficient  due  to  propeller) 
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ACt  = (ACt) J = f(CT)  (increase  of  local  skin  friction 

\Uo/  *-  p -l\uo/  coefficient  due  to  propeller) 


tan  a = dr/dx 

The  pressure  contribution  to  the  propeller-hull  interaction  is  approximated  by  the 
potential-flow  approximation  of  ACp.  The  propeller-induced  velocity  on  the  hull  can  be 
computed  from  the  propeller  field  point  program.31  To  determine  the  change  in  pressure  due 
to  the  presence  of  the  propeller,  the  Douglas-Ncumann  computer  program13  has  to  be 
employed  with  a nonzero  normal  velocity  at  each  quadrilateral  in  order  to  cancel  the  normal 
velocity  induced  by  the  propeller 
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3,Kerwin.  J.E.  and  R.  Leopold.  “A  Design  Theory  for  Subcavitating  Propellers,”  Transactions  of  the  Society 
of  Naval  Architects  and  Marine  Engineers,  Vol.  72,  pp.  294-315  (1964). 
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Here  n is  the  unit  outward  normal  to  the  surface,  vpr  is  the  propeller-induced  radial  velocity, 
and  u is  the  propeller-induced  axial  velocity.  Thus,  when  the  velocity  fields  of  the  hull  and 
the  propeller  are  superimposed,  the  normal  velocity  will  be  zero  along  the  body  contour.  The 
total  pressure  coefficient  in  the  presence  of  the  propeller  is 


(Cpl  = 
p P 


_^  = 1 - (1  + 
U,2 


sp 


U., 
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Here  Up  is  the  total  surface  potential  velocity  in  the  presence  of  the  propeller,  which  consists 
of  us  (the  perturbation  velocity  due  to  the  bare  hull).  usn  (perturbation  velocity  as  the  result 
of  cancelling  the  norma)  propeller-induced  velocity),  and  uip  = ua  cos  a + vpr  sin  a (the 
propeller-induced  tangential  velocity  on  the  hull).  The  original  pressure  coefficient  of  the 
bare  hull  is  denoted  by  (Cp >b-  Since  the  pressure-defect  component  of  thrust  deduction  is  of 
a potential-flow  nature,  the  computed  t due  to  the  interaction  of  propeller-induced  velocity 
on  the  hull.  Equation  (18),  is  therefore  equal  to  the  Beveridge1  formulation  of  the  effect  of 
hull-induced  velocity  on  the  propeller. 


PROPELLER  AND  STERN  BOUNDARY-LAYER  INTERACTION 


The  boundary-layer  velocity  profile  in  the  absence  of  the  propeller  is  called  the  nominal 
profile.  The  nominal  velocity  distribution  at  the  propeller  plane  of  the  ship  model  is  usually 
measured  by  a standard  wake  rake  before  the  actual  design  ot  the  lull-scale  propeller  is  under- 
taken. The  suction  of  the  propeller  accelerates  the  flow  in  the  vicinity  of  the  propeller.  The 
effective  velocity  profile  to  be  used  in  the  design  of  a wake-adapted  propeller  is  the  result  of 
the  interaction  of  the  propeller  and  the  thick  stern  boundary  layer.  The  radial  distribution  of 
effective  velocity  profile  is  sometimes  scaled  up  or  down  trom  the  measured  circumferential- 
mean  nominal  velocity  profile  by  a constant  tactor.  This  constant  factor  has  been  taken  as 
the  ratio  of  the  Taylor  wake  fraction  to  the  measured  volume-mean  nominal  wake  fraction. 
This  constant-factor  approach  can  be  considered  only  as  a rough  empirical  approximation. 

As  stated  earlier,  the  only  known  previous  etfort  to  address  this  problem  theoretically  was  by 
D.M.  Nelson*  who  developed  an  unpublished  computer  program  for  calculating  the  effective 
wake  from  the  measured  nominal  wake  and  static  pressure  distribution  across  the  boundary 
layer.  The  method  to  be  presented  below  was  developed  by  using  an  approach  different  from 
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that  of  Nelson  and  requires  only  the  measured  nominal  wake  distribution.  A serious  effort 
has  been  made  here  to  compare  the  theoretical  wake  predictions  with  wake  distribution 
measured  by  LDV  in  the  presence  of  a propeller. 

THEORETICAL  CONSIDERATIONS 

' The  experimental  evidence  shown  in  the  next  section  allows  one  to  conclude  that  the 
influence  of  the  propeller  on  the  axisymmetrical  stern  boundary  layer  is  limited  to  a region 
extending  two  propeller  diameters  upstream  of  the  propeller.  It  is  also  known,  from  theory, 
that  the  mean  circumferential  velocity  \q  upstream  of  the  propeller  is  identically  equal  to 
zero  with  (or  without)  the  propeller  in  operation.  The  following  assumptions  are  made  to 
derive  a theoretical  approximation  of  the  hydrodynamic  interaction  between  a propeller  and 
a fully  developed  thick  turbulent  boundary  layer  upstream  of  the  propeller.  ( 1 1 the  How  is 
axisyminetric  and  the  fluid  is  incompressible;  (2)  the  interaction  ol  propeller  and  stern 
boundary  layer  is  considered  to  be  inviscid  in  nature  and  thus  propeller-induced  viscous  losses 
and  turbulent  Reynolds  stresses  are  neglected;  (3)  the  conventional  boundary-layer  assumption 
of  9vr/9x  « 9ux/9r  is  considered  to  be  valid  for  the  nominal  boundary  layer  in  the  absence 
of  a propeller;  and  (4)  upstream  of  the  propeller,  no  energy  is  added  to  the  fluid  by  the 
propeller,  and  the  propeller-induced  velocity  field  upstream  ol  the  propeller  is  irrotational. 

The  vector  equation  of  steady  motion  for  an  inviscid  fluid  is  given  by  (see.  for  example. 
Thwaites32) 

V x oj  = - grad  H <21  > 

P 

_>  — ► — ► j —+ 

Here  V is  the  fluid  velocity,  u>  = V x V is  the  vorticity  vector,  and  H = P + -p(V*Vl  = 

P + i p(v2  + vy1  + ux2)  is  the  total  head.  For  cylindrical  polar  coordinates  (r.  0,  xl.  with 
V = <vr  v0,  ux).  the  radial  component  of  Equation  <21 ) may  be  written  in  the  form 

rvp  frrvfl)  _ | <jh  di b 

r2  9r  L*x  p dii/  dr 

where  the  total  pressure  head  H is  a function  only  of  the  stream  function  \j>  and 

= 9vf/9x  - 9ux'9r  is  the  0-component  of  the  vorticity  vector.  For  an  incompressible 
axisymmetrical  flow,  a stream  function  may  be  introduced  such  that 


^Thwaites.  H , “Incompressible  Aerodynamics.”  Oxford  University  Press  (1460),  Chapter  XI. 
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Since  upstream  of  the  propeller,  the  velocity  component  is  identically  zero  with  and  with- 
out the  propeller  in  operation.  Equation  (22)  reduces  to 

1 ™ = 1 (?!!*_  IM  (23) 

fj  di^  r \ 3r  3x  / 


upstream  of  the  propeller. 

Recall  that  the  boundary-layer  velocity  profile  in  the  absence  of  the  propeller  is  called 

the  nominal  profile.  The  suction  of  the  propeller  accelerates  the  flow  in  the  vicinity  of  the 

propeller  and  the  stream  surfaces  are  shifted  closer  to  the  body  surface.  As  shown  in  Figure  6, 

a typical  stream  surface  moves  inward  from  r to  rp  while  the  velocity  is  increased  from  the 

nominal  velocity  u to  u due  to  the  interaction  of  the  propeller  with  the  nominal  velocity 
x p 

profile.  TL  resultant  velocity  up  measured  i..  front  of  an  operating  propeller  will  be  called 
the  apparent  velocity  profile.  The  input  velocity  profile  to  be  used  in  the  design  of  a wake- 
adapted  propeller  will  be  called  the  effective  velocity  profile.  The  effective  velocity  profile  is 
the  apparent  velocity  profile  minus  the  irrotational-fiow,  propeller-induced  velocity  profile. 

*N 


Figure  6 - Illustration  of  Propeller/Stcrn  Boundary-Layer  Interaction 
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The  normal  boundary-layer  approximation  9vr/9x  « 9ux/9r  is  assumed  to  be  valid  for 
the  nominal  profile.  Therefore,  from  Equation  (23) 


I dH  1 /9ux  dvA  1 3ux 

p d^  = 7V"dT  ~ ar; -7  ~dT  wlthout  a propeller 


Applied  to  the  resultant  flow  with  the  propeller  in  operation.  Equation  (23)  results  in 
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Here  we  have  neglected  propeller  interaction  on  the  radial  velocity.  The  radial  velocity  with 
the  propeller  in  operation  is  equal  to  vf  + vpr,  where  vf  is  the  radial  velocity  without  propeller 
and  vpr  is  the  circumferential-average  propeller-induced  radial  velocity.  The  term  3vr/9x  used 
to  approximate  Equation  (24)  has  also  been  neglected  for  the  same  reason. 

Since  no  energy  is  added  to  the  upstream  fluid  by  the  propeller  and  since  we  assume  no 
viscous  losses  due  to  propeller-induction  effects,  the  total  pressure  head  within  the  same 
stream  annulus  remains  constant  with  and  without  the  propeller  in  operation,  i.e.,  upstream 
of  the  propeller 


i dH  = J_  /dH\ 

p d\p  p \d\p)p 

Furthermore,  we  have  assumed  that  the  propeller-induced  velocity  field  is  irrotational.  i.e.. 


V xVp 


where  Vp  is  the  total  propeller-induced  velocity,  and  vpr  and  ua  are  the  circumferential-average 
radial  and  axial  propeller-induced  velocities,  respectively.  With  these  two  conditions,  we 
obtain  a simple  formula  from  Equations  (24)  and  (25)  for  the  location  of  the  new  stream 
surface  rp 


1-  = 1 
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At  a given  location  of  x within  a given  stream  annulus  having  constant  mass  flow,  we  have  by 
definition 


d*  = ruxdr  = rP  VlrP 
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Thus,  Equation  (26)  may  be  rewritten  as 

uxdux  = Upd(Up  “V  (28) 

Equations  (27)  and  (28)  are  the  governing  equations  for  the  propeller  and  stern  boundary- 
layer  interaction.  (The  same  theory  is  also  valid  for  a windmilling  propeller  which  retards  the 
How  upstream  of  the  propeller.)  In  the  absence  of  a propeller,  the  nominal  velocity  profile 
(ux  versus  r)  and  the  propeller-induced  axial  velocity  profile  (ua  versus  r)  can  be  used  to 
obtain  the  new  location  of  the  stream  surface  rp  with  velocity  equal  to  up  via  Equations  (27) 
and  (28).  Equation  (27)  is  a statement  of  the  conservation  of  mass  in  an  incompressible  fluid: 
the  values  of  r and  rp  begin  at  r = rp  = rQ  (rQ  is  the  radius  of  the  body)  and  move  outward. 

Far  outside  the  boundary  layer  as  rp  -*•  dux  = 0 in  Equation  (28);  this  implies  d(up  - ua) 

= 0 and  u„  - u = constant  = U„  since  u„  is  zero  and  u„  has  to  be  equal  to  LL  as  r„  -*■  °°. 

P a U d P Up 

Thus,  for  rp  -*  «>,  the  effective  velocity  is  identically  equal  to  UQ  or  the  apparent  velocity 
u„  = LL  + u . These  results  are,  of  course,  valid  for  a propeller  operating  in  open  water  where 

p 0 3 

the  difference  in  the  nominal  velocity  dux  is  identically  equal  to  zero  since  ux  = U0  for  all  r. 
Thus,  the  effective  velocity  for  a propeller  operating  in  open  water  is  identically  equal  to  UQ 
even  though  the  apparent  velocity  up  is  equal  to  UQ  + ua(rp). 


NUMERICAL  SOLUTION 


The  nominal  velocity  ux  can  be  approximated  very  accurately  over  a small  increment  of 
radius  dr  by  a linear  function  of  r.  Although  the  velocity  at  the  wall  is  zero,  the  velocity 
profile  in  the  present  approximation  will  be  extrapolated  linearly  toward  the  wall,  resulting  in 
a nonzero  velocity  at  the  wall.  Since  uv  and  u_  can  be  approximated  locally  by  a linear 

X p 

function  of  r and  rp,  the  mass  flux  within  the  stream  tube  annulus  given  by  dr  = r;  + , - r;  and 
drp  = rpj  + j _ rpj  can  be  integrated  from  Equation  (27)  as 
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with  u(r)  = Uj  + [(r  - rt)/( r4  + , - r{)l  (uj  + , - u;)  and  r;<  r<  rj  + , 
to  obtain  the  finite  difference  form  as 


ri2+ , - fi2  ) (2  u„.  , + ux.)  - (ux.  - ux.) 
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= (r~  ~rD- ) (2u„  +uD.)-(uD  -Up.)  7 (29) 
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Equation  (28)  can  also  be  written  in  the  finite  difference  form 

(uxi  + , - % > = < uPi  + , + l‘Pi KuPi  + 1 ~ uPi  - % + , + V (30> 

Equations  (29)  and  (30)  are  simultaneous  equations  to  be  used  to  solve  for  rp  and  u 
A simple  computer  program  has  been  developed33  to  solve  them. 

PREDICTION  OF  PROPELLER  PERFORMANCE 

The  nominal  circumferential-mean  velocity  profile  at  the  propeller  plane  (x  = 0)  of  a 
ship  model  is  usually  determined  experimentally  from  pitot  tube  measurements  in  the  absence 
of  the  propeller.  The  effective  velocity  profile  is  the  profile  the  propeller  blades  experience 
in  producing  the  local  thrust  and  torque.  The  effective  profile  is  an  important  input  in  the 
design  of  the  modern  wake-adapted  propeller.  The  effective  profile  is  sometimes  scaled  up  or 
down  from  the  measured  nominal  profile  by  a constant  factor  which  is  taken  to  be  the  ratio 
(1  - wr)/(l  - w?),  where  wT  is  the  Taylor  wake  fraction  and  Wy  is  the  volumetric  mean 
nominal  wake  fraction  determined  from  a pitot-tube  survey.  An  accurate  effective  velocity 
profile  is  vital  for  the  successful  design  of  a wake-adapted  propeller  and  for  the  correct 
prediction  of  propeller  powering  and  cavitation  performance.  It  is  essential  to  develop  a 
reliable  and  sound  theoretical  procedure  to  calculate  the  effective  velocity  profile.  In  the 
present  investigation,  the  effective  velocity  profile  will  be  computed  by  a simple  iteration 
scheme  based  on  propeller/stern  boundary-layer  interaction  theory  and  a propeller  inverse 
program.34  A brief  description  is  outlined  below. 

The  nominal  wake  w(r)  is  defined  as 


- w(  r ) 


33Muang.  T.T.,  “User's  Manual  for  a Fortran  IV  Computer  Program  for  Calculating  Propeller,  Stem  Boundary- 
Layer  Interaction  on  Axisymmetrical  Bodies.”  DTNSRDC  Report  SPD-737-02  (1976). 

34Cummings,  D.F..  “Numerical  Prediction  of  Propeller  Characteristics.”  Journal  ot  Ship  Research.  Vol.  17. 

Pt  3.  pp.  12-18  (1973). 
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The  effective  wake  experienced  by  a propeller  blade  section  is  defined  by 
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In  current  design  practice,  the  effective  wake  distribution  is  sometimes  estimated  by  the 
following  approximation 
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The  Taylor  wake  fraction  wT  is  determined  experimentally  from  self-propulsion  and 
open-water  propeller  tests  in  a towing  tank;  w^  is  the  volumetric  mean  wake  derived  from 
wake  surveys  conducted  in  the  absence  of  a propeller.  Because  of  the  difference  in  the 
propeller  loading  and  ambient  turbulence  (which  affects  the  propeller  section  drag  coefficient), 
the  propeller  characteristics  in  the  behind  condition  can  be  quite  different  from  characteristics 
in  open  water.  Furthermore,  the  propeller  interaction  on  the  boundary  layer  cannot  be 
represented  by  a constant  correction  factor.  In  order  to  improve  the  prediction  of  the 
effective  velocity  profile  from  the  measured  velocity  profile  at  the  propeller  plane,  we  use  the 
following  simple  iteration  scheme  based  on  the  previously  discussed  propeller/stern  boundary- 
layer  iteraction  program  and  the  propeller  inverse  program  developed  by  Cummings:34 

1.  Use  the  measured  nominal  velocity  profile  to  estimate  an  effective  velocity  profile. 
One  may  start  with  the  guess  ug(r)  =1.1  ux(r). 

2.  Use  the  propeller  inverse  program  developed  by  Cummings34  to  calculate  the  propeller 
nondimensional  circulation  distribution  G(rp)*  from  the  estimated  effective  velocity  profile 
ue(rp)  and  the  propeller  geometry. 

3.  Use  the  propeller  field-point  velocity  program31  to  calculate  the  average  propeller- 


induced  axial  velocity  ua. 


“The  inverse  program  has  to  be  calibrated  with  the  conventional  design  program31,35'36  for  the  given 
effective  profile  for  which  the  particular  propeller  is  designed. 

35Cox.  G.G.  and  W.B.  Morgan.  “Use  of  Theory  in  Propeller  Design,”  Marine  Technology.  Vol.  2,  No.  4, 
pp.  319-320  (1072). 

36 Morgan,  W.B  et  al„  “Propeller  Lifting-Surface  Corrections.”  Transactions,  the  Society  of  Naval  Architects 
and  Marine  Engineers,  Vol.  76.  pp.  300-347  (1068). 
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4.  Use  the  propeller/stern  boundary-layer  interaction  program  to  compute  up( rp ) from 
the  measured  value  of  uv(r)  and  computed  value  of  u (r  ). 

A dp 

5.  The  new  effective  profile  becomes  ue(rp)  = up( rp ) - ua( rp ). 

6.  Repeat  Steps  2-4  until  uJr_)  and  G(r_)  converge. 

c p p 

The  final  computed  nondimensional  circulation  distributions  G(r)  for  the  three  afterbodies 
are  shown  in  Figure  7.  It  has  been  found  that  the  values  of  ue/UQ  and  G(r)  have  essentially  con- 
verged to  their  final  values  after  three  iterations.  The  computations  also  show  that  the  ratio  of  the 
effective  to  nominal  velocity  distribution  t*e(rp^ux^rp * = up(rp)/ux(rp) ~ ’ " ua<rp^ux<rp* 
does  not  take  a constant  value.  In  fact,  the  ratio  of  u.(r_)/uv(r_)  has  a larger  value  near  the 

C p A p 

hub  of  the  propeller  than  near  the  tip. 


COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  RESULTS 

MEASURED  AND  COMPUTED  PRESSURE 
AND  SHEAR  STRESS  DISTRIBUTIONS 


The  measured  and  computed  values  of  pressure  coefficient  are  compared  in  Figure  8 for 
Afterbodies  1,  2,  and  3 at  a length  Reynolds  number  Rn  = 5.9  x I06.  The  virtual  origin  of 
turbulence  was  fixed  at  X/L  = 0.015  (Figure  4b).  The  agreement  was  very  good  for  Afterbody 
1,  fairly  good  for  Afterbody  2.  and  relatively  poor  for  Afterbody  3.  Recall  that  Afterbody  1 
did  not  separate.  Afterbody  2 approached  separation,  and  Afterbody  3 separated  at  X/L  = 0.42. 
Because  of  the  severe  adverse  pressure  gradients  on  Afterbodies  2 and  3,  the  boundary  layer 
became  very  thick  and  the  accuracy  of  the  predictions  was  degraded.  The  results  suggest  that 
the  present  method  is  accurate  for  sterns  without  severe  adverse  pressure  gradients  but  less 
accurate  for  sterns  with  progressively  more  severe  pressure  gradients.  Thus,  further  improve- 
ment of  the  present  computation  method  is  necessary. 

Figure  4 shows  the  measured  and  computed  difference  between  the  viscous  How  and 
potential  flow  pressure  (Cpv(-  -Cpp)  for  the  three  afterbodies.  The  shaded  areas  indicate  the 
additional  drag  caused  by  separation  on  Afterbody  3 and  the  additional  measured  stern 
pressure  drag  over  the  prediction  on  Afterbody  2. 

The  measured  and  computed  local  shear  stress  distributions  C'T  are  compared  for  After- 
body I in  Figure  10.  for  Afterbody  2 in  Figure  1 1.  and  for  Afterbody  3 in  Figure  12.  Two 
body-length  Reynolds  numbers  (Rn  = U0L/ri  of  5.4  x I0h  and  8.8  x lO6  were  used  in  the 
experiments.  As  can  be  seen  from  these  figures,  the  agreement  was  good  for  Afterbodies  1 
and  2.  However,  for  Afterbody  3,  agreement  between  measured  and  computed  shear  stresses 
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Figure  9 — Measured  and  Computed  Viscous  Effect  on  the  Stern 
Pressure  Distributions  of  the  Three  Afterbodies 


was  good  only  up  to  the  point  of  separation  (X/L  = 0.92);  after  separation,  the  present 
computational  method  breaks  down  and  the  Preston  tube  measurement  is  also  questionable. 
Therefore,  Figures  12a  and  12b  do  not  include  results  after  separation.  Qualitatively,  the 
shear  stress  takes  small  negative  values  in  the  separation  bubble  and  becomes  positive  again 
after  flow  reattachment  at  about  X/L  = 0.97.  It  should  be  pointed  out  that  the  computed 
values  of  CT  for  Afterbodies  1 and  2 were  higher  than  the  measured  values  of  C'r  at 
X/L  > 0.95.  This  discrepancy  is  similar  to  that  for  the  measured  and  computed  pressure 
distributions.  The  disagreement  between  measured  and  computed  values  of  C'T  at 
X/L  > 0.95  was  more  noticeable  for  Afterbody  2 than  Afterbody  1. 
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Figure  12  — Measured  and  Computed  Skin  Friction  Distribution 
on  \fterbodv  with  and  without  Propeller  in  Operation 


*The  effect  of  the  propeller  on  the  stern  pressure  distribution  was  computed  by  the 
potential-flow  approximation  method  described  earlier.  As  can  be  seen  from  Equation  (20). 
the  total  surface  perturbation  potential  velocity  in  the  presence  of  the  propeller  is  the  sum  of 
the  bare-hull  perturbation  velocity  (including  the  displacement  effect),  the  perturbation 
velocity  as  the  result  of  cancelling  the  propeller-induced  velocity  normal  to  the  hull,  and  the 
propeller-induced  velocity  tangential  to  the  hull.  The  boundary-layer  calculation  can  be  made 
by  using  the  computed  total  surface  perturbation  potential  velocity  in  the  presence  of  a 
propeller. 

In  connection  with  the  measured  and  computed  values  of  C_  with  the  propeller  in 
operation  (Figures  10-  12),  note  that  as  expected,  the  propeller  accelerated  the  How  at  the 
stern,  resulting  in  an  increase  of  shear  stress.  However,  this  increase  of  Cr  was  limited  to  the 
region  X/L  > 0.90.  No  effect  was  noted  at  distances  larger  than  two  propeller  diameters 
(2Dp)  upstream  from  the  propeller.  As  shown  in  Figures  12a  and  12b,  the  suction  of  the 
propeller  did  not  change  the  point  of  boundary-layer  separation  on  Afterbody  3.  The 
distance  between  the  propeller  plane  and  the  point  of  separation  was  1.3  propeller  diameters. 
The  propeller-induced  velocity  at  1.3  I)p  upstream  of  the  propeller  was  not  strong  enough  to 
alter  the  characteristics  of  the  separation. 

Figures  13-15  show  the  measured  and  computed  pressure  distributions  with  and  without 

the  propeller  in  operation.  Again,  the  effect  of  the  propeller  on  pressure  distribution  was 

evident  up  to  a distance  of  2Dp  upstream  of  the  propeller.  As  can  be  seen  from  Figures  13 

and  14,  the  measured  values  of  (Cn)  , (Cn),  , and  AC'  = (Cn),  -<Cn)  were  in  good  agreement 

P p P b P P b P p 

with  the  computed  values  for  Afterbodies  1 and  2.  In  the  case  of  Afterbody  3 (Figure  15), 
the  measured  values  of  ACp  were  smaller  than  the  computed  values  of  ACp  aft  of  the 
separation  point  X/L  = 0.92.  These  smaller  measured  values  of  ACp  may  have  been  caused 
by  the  cushion  effect  of  the  separation  bubble.  Upstream  of  the  separation  point,  however, 
the  measured  values  were  larger  than  the  computed  values  of  ACp.  reflecting  the  possible 
contraction  of  the  separation  streamline  due  to  the  propeller. 

MEASURED  AND  COMPUTED  THRUST  DEDUCTION 

As  shown  in  Figures  13  and  14.  the  measured  values  of  <C’n)  . <C'),  , and 

P p P b 

A(  „ = (('_),  — < C’  > were  in  good  agreement  with  the  computed  values  for  Afterbodies  1 and 
P P b P p 

2.  The  agreement  for  AC'p  was  less  satisfactory  for  Afterbody  3 near  the  separation  point  of 
X/L  = 0.92  (Figure  15).  These  figures  also  show  the  values  of  tp  calculated  from  the 
integration  of  the  measured  values  ot  A('p  by  using  Equation  ( IX). 
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The  good  agreement  between  the  measured  and  computed  values  of  ACp  shown  in 
Table  4 for  Afterbodies  1,  2,  and  3 indicates  that  the  pressure  component  of  thrust  deduction 
can  be  well  predicted  by  the  potential-flow  propeller-hull  interaction  approximation.  However, 
if  there  is  flow  separation  at  the  stem,  the  potential-flow  interaction  approximation  for  tp 
should  be  used  with  caution.  Although  the  computed  value  of  t for  Afterbody  3 agreed 
well  with  the  measured  value  of  tp  obtained  from  the  interaction  of  the  measured  ACp,  the 
latter  did  not  agree  very  well  with  the  computed  ACp  in  the  vicinity  of  the  separation  point 
(Figure  15). 

The  CS  boundary-layer  computation17  can  be  performed  with  the  new  pressure  distribu- 
tion (Cp)p.  The  measured  and  computed  local  shear  stress  distributions  and  the  value  of  tF 
computed  by  integration  of  measured  values  of  ACTo  (by  using  Equation  (19))  over  the  three 
afterbodies  have  already  been  shown  in  Figures  10-12.  The  measured  value  of  t and  tf;  for 
the  three  afterbodies  at  Rn  = 5.9  x 106  are  summarized  in  Table  4. 

TABLE  4 - COMPARISON  OF  COMPUTED  AND 
MEASURED  THRUST  DEDUCTION 


Configuration 

Computed* 

lP 

Measured 

lP 

Measured 

lF 

Measured 

t;MP+tF 

Measured 

Afterbody  1 , J = 1 .25 

0.068 

0.07 

0.0024 

0.0724 

3.4% 

Afterbody  2,  J = 1 .25 

0.129 

0.143 

0.0023 

0.1453 

1.6% 

Afterbody  2.  J = 1.07 

0.126 

0.140 

0.0026 

0.1426 

1.9% 

Afterbody  3,  J = 1.25 

0.106 

0.109 

0 

0.109 

0 

Afterbody  3,  J = 1.07 

0.103 

0.103 

0 

0.103 

0 

*The  computed  values  of  tp  were  obtained  by  integration  of  the  computed  values  of  ACp. 


Agreement  was  satisfactory  between  the  computed  and  measured  pressure  components 
of  thrust  deduction  tp.  The  measured  ratio  tf  /tp  was  less  than  5 percent  for  the  three  after- 
bodies at  J = I 25  and  J = 1.07.  Therefore,  the  value  of  thrust  deduction  t can  be  satisfactorily 
estimated  from  its  pressure  component  tp  and  its  frictional  component  tr  can  be  neglected. 

MEASURED  AND  COMPUTED  AXIAL  VELOCITY  PROFILES 

The  mean  axial  velocity  profiles  in  the  boundary  layers  of  the  three  afterbodies  with  and 
without  a propeller  in  operation  were  measured  by  a laser  doppler  velocimetcr.  The  boundary- 
layer  profiles  were  also  computed  bv  the  Douglas  CS  differential  boundary-layer  computer 
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program.13,17'30  The  pressure  distributions  calculated  from  the  present  potential-flow 
boundary-layer  iteration  scheme  were  used  to  compute  the  protiles  tor  the  bare  hull.  On  the 
other  hand,  the  pressure  distributions  calculated  from  the  potential-flow  propeller-hull 
interaction  approximation.  Equation  (20),  were  used  to  compute  the  protilcs  with  propeller 
in  operation.  The  measured  and  computed  mean  axial  velocity  profiles  are  compared  in 
Figures  16  to  19.  The  present  approximations  for  the  profiles  with  a propeller  in  operation 
were  computed  from  the  inviscid  propeller/stern  boundary-layer  interaction  approximation. 

Figures  16a-  16e  respectively  show  the  measured  and  computed  axial  velocity  protiles 
u /U  for  Afterbody  1 at  X/L  = 0.845,  0.915.  0.954,  0.964,  and  0.977.  The  agreement  was 
reasonably  good  except  at  X/L  = 0.977  where  the  measured  velocities  were  noticeably  smaller 
than  the  computed  velocities  near  the  body  in  the  absence  of  a propeller.  The  generally  good 
agreement  between  measured  and  computed  velocity  profiles  (except  very  close  to  the  body) 
on  Afterbody  1 is  consistent  with  the  good  agreement  between  the  measured  and  computed 
shear  stress  (Figure  10)  and  pressure  distributions  (Figure  8). 

Corresponding  data  for  Afterbody  2 at  X/L  = 0.935,  0.964.  and  0.677  are  shown  in 
Figure  17.  Agreement  was  good  at  X/L  = 0.935  but  the  measured  velocities  were  progressively 
slower  than  the  computed  velocities  as  the  end  of  the  stern  was  approached.  This  discrepancy 
near  the  stern  is  consistent  with  that  found  between  the  measured  and  the  computed  shear 
stresses  (Figure  I I ) and  pressure  distribution  (Figure  8)  on  Afterbody  2.  The  present 
potential-flow  boundary-layer  iteration  scheme  is  not  adequate  for  handling  the  How  on  a full 
afterbody  with  a small  value  of  LA/D  (2.24)  and  with  a strong  adverse  pressure  gradient. 
However,  the  iteration  scheme  does  provide  a very  good  approximation  of  the  stern  flow  for 
Afterbody  1 which  has  a mild  adverse  pressure  gradient  (Figure  2i  and  a moderate  value  of 
L4/D  (=  4.3). 

Figure  18  shows  only  the  measured  axial  velocity  protiles  at  X 1.  = 0.963  and  0.677  tor 
Afterbody  2 with  the  propeller  moved  one-quarter  propeller  diameter  ( D(1  4)  att  of  its 
original  location  at  X/L  = 0.983.  Comparing  Figures  17b  and  17c  with  1 8a  and  Figures  Pd 
and  1 7e  with  18b.  one  finds  as  expected  that  the  velocity  augmentation  due  to  propeller 
suction  is  reduced  by  moving  the  propeller  further  downstream. 

Figures  19a-  19e  respectively  show  the  measured  and  computed  velocity  profiles  for 
Afterbody  3 at  X/L  = 0.795  and  0.6)5  and  only  measured  values  at  X L = 0.934.  0.954.  and 
0.677.  The  agreement  was  excellent  at  X/L  = 0.76.5  (Figure  1 6a)  but  only  fair  at  X L = 0.615 
(Figure  19b)  immediately  upstream  of  the  separation  point  (X  L ~ 0.62).  It  is  also  important 
to  note  tiiat  the  measured  axial  velocity  approached  the  tree-stream  velocity. 
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Figure  16  — Measured  and  Computed  Axial  Velocity  Profiles 
on  Afterbody  1 with  and  without  Propeller  in  Operation 
at  Different  X/L  Ratios  and  with  R„  = 5.9  x 106 
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Figure  16c  - At  X/L  = 0.954 
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Figure  17b  - At  X/L  = 0.964,  B-L  Calculation 
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Figure  18b  - At  X/L  = 0.977 

Figure  18  - Measured  Axial  Velocity  Profiles  on  Afterbody  2 
with  Propeller  Moved  Dp/4  Aft  of  Original  Location 
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Figure  19e  - At  X/L  = 0.977 

On  the  other  hand,  within  the  boundary-layer  approximation,  the  computed  axial  velocity 
approached  the  local  potential  axial  velocity  on  the  body  at  the  edge  of  the  boundary  layer. 
The  separated  How  reattached  at  X/L  = 0.97.  Since  the  CS  boundary-layer  computer  program 
terminates  computation  when  How  separation  is  encountered,  no  computed  results  are 
presented  at  X/L  = 0.934,  0.954,  and  0.977.  The  effect  of  the  propeller  on  the  measured 
profiles  at  X/L  = 0.934  and  0.954  (Figures  19c  and  19d)  was  relatively  small,  but  it  became 
significant  at  X/L  = 0.977  (Figure  19e)  where  the  flow  reattached  after  separation. 

The  fact  that  the  effect  of  the  propeller  on  the  measured  velocity  profiles  of  the  three 
afterbodies  was  limited  to  distances  of  two  propeller  diameters  upstream  of  the  propeller 
(Figures  18  and  19)  was  consistent  with  the  measured  effect  of  the  propeller  on  shear  stress 
(Figures  10- 12)  and  pressure  (Figures  1 3 - 1 5)  distributions. 

Let  us  briefly  review  data  already  presented  for  the  measured  and  computed  axial  velocity 
profiles  u(1(r)  with  the  propeller  in  operation.  Fxamine  Figure  16c  (X/L  = 0.954)  and 
Figure  I6e  (X/L  = 0.977)  for  Afterbody  1.  Figures  17b  and  17c  (X/L  = 0.9(i4l  and  Figures 
I 7d  and  I 7e  (X/L  = 0.977)  for  Afterbody  2.  and  Figure  19e  (X/L  = 0.977)  for  Afterbody  3. 
Agreement  between  the  measured  and  computed  values  of  u^  was  good  for  all  of  the  cases. 

I(  is  important  to  note  that  the  propeller  plane  was  at  X/L  = 0.983  in  all  cases,  lhc  measure- 
ment stations  were  located  at  small  distances  from  the  propeller.  The  propeller  field-pr  int 
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velocity  program31  was  used  to  calculate  the  circumferential-mean,  propeller-induced  axial 
velocities  at  the  positions  where  the  nominal  and  effective  velocity  profiles  were  measured. 

The  change  in  the  radial  distribution  of  circulation  for  the  propeller  in  the  wakes  of  the  three 
different  afterbodies,  needed  for  the  propeller  field-point  velocity  program.31  has  already  been 
discussed. 

The  propeller  field-point  velocity  program31  has  two  options:  one  uses  lifting-line  theory 
md  the  other  uses  lifting-surface  theory.  The  mean  propeller-induced  axial  velocities  ua 
calculated  by  the  two  options  differ  somewhat  near  the  propeller  but  the  difference 
diminishes  as  distance  from  the  propeller  is  increased.  Figures  20-22  present  the  computed 
profiles  ua/U()  at  various  measurement  stations.  The  difference  in  axial  velocity  with  and 
without  a propeller  is  defined  as 


U 


o 


II 


u 


The  measured  and  calculated  values  of  Aux/U0  at  the  various  stations  of  the  three  after- 
bodies are  also  shown  in  Figures  20-  22.  As  indicated  in  Figure  20a  < X I 0 >,>4.  0.5H  Dp 
upstream  of  the  propeller),  the  measured  and  calculated  values  of  Aux  l 0 were  in  good 
agreement  and  this  was  so  even  for  the  values  of  Aux/U0  calculated  by  the  boundary-layer 
method,  with  the  wall  pressure  distribution  modified  to  account  tor  the  suction  of  the 
propeller.  However,  at  X/L  =0.977  (0.12  Dp>,  the  values  of  Aux/UQ  determined  by  the 
boundary-layer  calculation  did  not  agree  well  with  the  measured  data  (Figures  20b.  21c.  and 
2ld>.  The  present  approximation  with  the  lifting-surface  option  was  in  good  agreement  with 
all  the  measured  data. 

Thus,  it  may  be  concluded  that  the  present  inviscid  approximation  provides  a good 
representation  of  the  complex  hydrodynamic  interaction  between  the  propeller  i.id  the  stern 
boundary  layer.  The  good  agreement  between  the  measured  and  computed  values  of  Aux/U0 
at  small  distances  from  the  propeller  suggests  that  the  present  approximation  can  be  used  with 
confidence  to  calculate  the  effective  velocity  profile  for  the  propeller  from  the  measured 
nominal  velocity  profile  and  a computed  average  propeller-induced  axial  velocity  profile. 


Figure  20a  - At  X/L  = 0.954 
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Figure  20b  - At  X/L  = 0.977 

Figure  20  - Measured  and  Computed  Axial  Velocity  Increase  on  Afterbody 
Due  to  Propeller  Suction  at  Different  X/L  Ratios 


Figure  21  — Measured  and  Computed  Axial  Velocity  Increase  on  Afterbody  2 
Due  to  Propeller  Suction  at  Different  X/L  Ratios 


Figure  21a  - At  X/L  = 0.964,  Propeller  J = 1.25 


Figure  21b  - At  X/L  = 0.964,  Propeller  J = 1.07 
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CONCLUSIONS 

A comprehensive  set  of  experimental  data  has  been  presented  for  three  axisymmetrical 
afterbodies,  two  with  attached  flows  and  one  with  separated  How.  In  addition  to  enabling 
analytical  techniques  to  be  evaluated,  these  data  also  serve  to  give  insight  into  the  details  of 
the  flow  in  the  stern  region.  The  following  conclusions  were  drawn: 

1.  The  potential  flow  boundary-layer  interaction  program  computed  accurate  values  of 
pressure  and  shear  stress  on  Afterbody  1 (which  has  a fine  stern),  gave  less  satisfactory 
predictions  for  Afterbody  2 (which  has  a fuller  stern)  and  were  in  only  fair  agreement  for 
Afterbody  3 (which  separates  at  X/L  = 0.92). 

2.  For  all  sterns,  the  agreement  between  computed  and  measured  velocity  profiles  was 
excellent  up  to  about  X/L  = 0.9.  Thereafter,  the  agreement  progressively  worsened  as  the 
tail  was  approached,  particularly  for  Afterbody  3.  These  results  suggest  that  the  modeling  of 
the  wake  should  be  modified  to  handle  cases  of  separated  flow  and  that  the  boundary-layer 
equations  should  be  modified  to  properly  model  the  thick  stern  boundary  layer. 

3.  In  most  cases,  the  propeller  stern  boundary-layer  interaction  program  predicted 
reasonably  accurate  effective  velocity  profiles.  The  difference  between  the  computed  effective 
and  nominal  velocity  profiles  was  largest  at  the  propeller  hub  and  became  smaller  at  the 
propeller  tip.  This  is  in  contrast  to  results  obtained  by  the  conventional  method  wherein  the 
effective  wake  is  scaled  up  or  down  from  the  nominal  wake  by  a constant  factor,  a procedure 
which  usually  yields  the  largest  difference  between  the  effective  and  nominal  velocity  profiles 
at  the  propeller  tip.  Since  the  present  experimental  results  compared  well  with  the  present 
theory,  it  is  concluded  that  the  constant-factor  method  is  not  satisfactory  and  should  no 
longer  be  used  for  bodies  of  revolution.  The  present  theory  takes  as  input  the  measured 
nominal  velocity  profile  without  the  propeller  in  operation  and  uses  standard  propeller 
computer  codes  for  computing  propeller  blade  loading  and  induced  field-point  velocities. 

4.  The  use  of  lifting-surface  theory  to  predict  the  induced  field-point  velocities  usually 
gives  more  accurate  effective  velocity  profiles  than  provided  by  lifting-line  theory. 

5.  When  propeller-induced  pressure  distribution  is  used  to  calculate  velocity  profiles 
near  the  propeller  by  the  CS  boundary-layer  method,  predictions  usually  deviate  substantially 
from  the  measured  velocity  profiles  in  the  presence  of  a propeller. 

6.  The  experimental  results  show  that  the  potential-flow  propeller-hull  interaction 
methods  predict  thrust  deduction  and  the  propeller-induced  pressure  distribution  very  well. 
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APPENDIX 

OFF-AXIS  DUAL-BEAM  BACKSCATTER 
LASER  DOPPLER  VELOCIMETER 


A focal  length  of  1.5  m was  selected  for  the  LDV  optics  so  that  the  entire  optical  bench 
could  be  located  in  the  quiescence  region  of  the  anechoic  chamber.  The  selection  of  a beam 
separation  of  99  mm  with  a total  beam  angle  of  3.722  deg  was  based  on  the  desu  d size  of 
probe  volume. 

As  shown  in  Figure  23,  the  collimated  beam  emitted  by  a multicolor  argon-ion  laser 
(Coherent  Radiation  Model  CR-3)  is  tuned  to  emit  green  (X  = 5145  A°t  light  and  is  polarized 
parallel  to  the  x-axis  by  a polarization  rotator.  Such  polarization  is  necessary  to  ensure 
maximum  scatter  intensity.37  A beam  splitter  separates  the  beam  into  two  equal-intensity 
parallel  beams  with  a separation  of  50  mm.  The  beams  are  separated  further  by  prisms  to  a 
final  separation  of  99  mm.  A 1 5.24-cm-diameter  lens  is  used  to  focus  the  two  parallel  beams 
at  probe  volume.  The  local  diameter  for  a diffraction-limited  system  is  given  38-39  by 


Here  b is  the  radius  of  the  focused  beam,  b is  the  radius  of  the  unfocused  beam,  X is  the 
wave  length  of  the  laser,  and  FL  is  the  focal  length  of  the  focusing  lens.  Since  each  beam 
can  be  assumed  to  have  a Gaussian  intensity  distribution,  the  probe  volume  will  be  ellipsoidal 
in  shape  and  have  the  following  dimensions: 


AX 


0.491  mm 


= 0.491  mm 
= 1 5.09  mm 


- 7 Born.  M.  and  E.  Wolf,  “Principles  of  Optics."  Pergarnon  Press,  New  York  (1965).  p.  652. 

38Coethert.  W.H.  and  D.B.  Brayton.  “New  Velocity  Measuring  Technique  Using  Dual  Scatter  Laser  Doppler 
Shift."  Arnold  Engineering  Development  Center  AFDC-TR-205  (1970). 

39Vanta.  W.J.,  “Turbulence  Measurements  with  a Laser  Doppler  Velocimeter Naval  Ordnance  Laboratory 
NOl.TR-73-94  (1973). 
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Here  AX  is  in  the  direction  of  the  mean  flow,  AY  is  perpendicular  to  the  mean  flow,  and 
AZ  is  in  the  bisector  of  the  beam  angle.  The  parameters  for  the  laser  and  optics  are  as 
follows: 

X = 5145  A° 

2 b = 2.0  mm  at  l/'e2  with  beam  expander 
FL  = 1500  mm 
0 = 3.722  deg 

The  large  size  of  AZ  is  primarily  due  to  the  beam  angle  6 which  is  the  ratio  of  lens 
diameter  to  its  focal  length.  The  present  choice  of  lens  diameter  (15.24  cm)  was  limited 
mainly  by  expense  since  the  cost  of  the  lens  increases  drastically  with  increasing  diameter. 
The  back-scattered  radiation  was  collected  vertically  off  axis  (=  28  deg)  by  a 15.24-cm- 
diameter  lens  which  focuses  the  scattered  light  onto  a photomultiplier  tube.  This  off-axis 
dual-beam  back  scatter  mode  of  operation  was  chosen  to  reduce  the  background  noise  from 
the  radiation  of  the  two  long  transmitting  beams;  it  can  also  be  used  to  reduce  the  effective 
length  of  the  long  axis  of  probe  volume  AZ.  An  effective  AZ  of  5 mm  can  be  achieved  by 
carefully  focusing  the  receiving  optics  on  the  center  portion  of  the  volume. 


BACKSCATTER 

MODULE 


In  view  of  the  coherence  characteristics  of  the  laser  light,  the  two  beams  will  interface 

0 

constructively  and  destructively  to  establish  a set  of  closely  spaced,  planar  interference  fringes 
at  the  probe  volume.  The  peak-to-peak  fringe  spacing  dQ  is 


d 


o 


= 7.9  pm 


The  number  of  fringes  contained  in  a probe  volume  can  be  determined  by  dividing  the  peak- 
to-peak  fringe  spacing  into  the  focal  diameter  AX  as  follows: 


. AX  4bo  le\  4 2FL  ie\  4 

N = r‘T,a"u)'»  !r,anbri 


beam  separation 
2b 


where  beam  separation  = 2FL  tan  (6/2)  = 99  mm  and  2b  = 2 mm.  Thus,  the  total  number 
of  fringes  in  the  probe  volume  is  N = 63,  which  is  sufficient  for  the  present  pc  nose  (N  must 
be  greater  than  10). 

As  a particle  traverses  the  probe  volume  at  a velocity  Ux,  the  interference  fringes  are  cut 
at  a rate  of 


f 


„ 9 

U sin  -z 


f = 


or  Ux 


T 


where  T is  the  period  of  Doppler  frequency  f.  The  unique  feature  of  this  dual-beam  scatter 
technique  is  that  it  is  possible  to  place  the  receiving  optics  in  any  position  since  the  view 
direction  does  not  affect  the  measurement.  It  should  be  noted  that  forward  scatter  has  a 
higher  signal-to-noise  ratio:  however,  the  mechanical  difficulty  associated  with  placing  the 
transmitting  and  receiving  optics  3 m apart  on  an  optical  bench  prevented  the  use  of  tcrward 
scatter  techniques  for  this  experiment.  The  off-axis  dual-beam  backscatter  technique  was 
used  as  a tradeoff  between  mechanical  difficulty  and  optical  difficulty. 

In  order  to  produce  satisfactory  backscatter  radiation,  it  was  necessary  to  seed  the 
boundary  layer  with  a fine  oil  mist;  the  mist  was  generated  by  a standard  atomizer  with  a 
Laskin  nozzle and  oil  was  emitted  from  an  orifice  it  the  model  bow.  The  size  ot  the 
particles  generated  varied  between  1 to  10  pin.  peaking  at  about  3 pm;  this  was  very  satisfac- 
tory compared  to  the  tringe  spacing  of  7.*)  pm. 

The  backscatter  radiation  of  these  particles  cutting  through  the  fringes  was  detected  by  a 
photomultiplier  operated  at  a 1400-Y  potential.  The  signals  from  the  photomultiplier  were 
filtered  and  amplified,  and  were  then  processed  h'_  i laser  Doppler  digital  data  processor 
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(Arnold  Research  Organization,  Inc.  Model  8)  which  has  been  described  in  detail  by  Kalb.40 
The  raw  data  were  recorded  on  a Hewlett-Packard  Model  505 5 A Digital  Printer,  which  allowed 
a maximum  data  rate  of  10  points  per  second. 

For  completeness,  a general  operating  procedure  will  be  given  for  the  counter-type  digital 
data  processor  (DDP).  Before  the  Doppler  frequency  or  “frequency-burst”  signal  can  be 
processed,  several  tasks  must  be  performed.  The  “pedestal”  voltage  and  other  gross  low- 
frequency  baseline  shifts  must  be  filtered  out  along  with  noise  above  and  below-  the  frequency 
range  of  interest.  The  pedestal  voltage  arises  because  of  unequal  beam  intensities,  or  polariza- 
tion of  the  two  beams,  or  because  the  particle  does  not  pass  through  the  geometrical  center 
of  the  scattering  volume.  This  pedestal  voltage  is  removed  by  a band-pass  filter  applicable  to 
the  signal  frequency  of  interest.  The  nine  bands  of  the  DDP  each  cover  2 1/2  octaves: 
together  they  provide  wide  frequency  coverage  from  15  kHz  to  50  MKz.  An  oscilloscope- 
triggered  gate  is  used  to  synchronize  the  processor  sampling  interval  with  the  burst  signal 
event  displayed  on  the  oscilloscope. 

The  DDP  employs  high-speed,  emitter-coupled  integrated  circuit  flip-flops  to  convert  the 
first  eight  pulses  of  Doppler  frequency  data  into  a square  waveform  with  the  same  frequency 
as  the  original  burst  signal.  The  processor  then  digitizes  the  resultant  time  interval  via  period 
counter  techniques.  The  DDP  then  uses  high-speed  digital  computing  circuitry  to  test  the 
data  pulse  train  for  periodicity.  It  compares  the  time  interval  of  both  four  and  five  data 
pulses  to  the  time  interval  of  eight  data  pulses.  The  two-stage  time  interval  comparison 
enhances  the  accuracy  of  the  Doppler  readings,  especially  when  sampling  low  signal-to-noise 
data.  The  use  of  both  4/8  and  5/8  time  interval  comparators  significantly  reduces  alias  or 
false  readings  caused  by  pulse  dropouts  and  pulse  additions.  The  DDP  has  a switchable  error 
band  available  with  0.78,  1.5,  and  3.0  percent  for  both  4 t and  5/8  time  interval  comparators. 
If  the  errors  of  the  two  comparators  are  within  the  error  band  chosen,  the  data  are  accepted 
for  recording  in  the  form  of  natural  binary  for  a minicomputer  or  magnetic  tape  recorder  or 
in  the  form  of  binary  code  decimal  for  conventional  digital  printers. 

After  validation,  data  are  accepted  by  the  data  acquisition  system  which,  in  turn,  releases 
the  recycle  sequence  switch  and  issues  a zero  reset  pulse  to  all  counters.  At  reset  pulse 
termination,  the  processor  must  further  await  the  reception  ol  a +A  gate  signal  trom  the 
oscilloscope  before  a new  sample  interval  may  begin.  The  recycle  sequence  is  altered  when 
the  errors  of  one  or  both  comparators  exceed  the  error  band  chosen.  \ “data  reject  signal 


4(,Kalb.  H.T.  et  al.,  “laisei  Velocimetry  Data  Processing,’'  Arnold  Engineering  Development  Center 
Ati>C-TR-73  116  (!‘>73). 
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will  then  be  issued  by  the  comparator.  All  counters  are  reset  to  zero,  the  invalid  data  will  be 
inhibited,  and  the  data  will  not  be  transferred  to  the  data  acquisition  system.  A new  data 
processing  cycle  will  be  initiated. 

The  valid  data  from  DDP  were  printed  by  a digital  printer  (Hewlett-Packard  Model  5055A) 
at  a rate  of  about  10  points  per  second.  At  each  profile  location  of  the  boundary  layer.  200 
to  400  Doppler  data  points  were  taken  and  used  as  input  to  the  CDC  6700  digital  computer. 
The  instantaneous  velocity  uxi  was  calculated  from  the  instantaneous  Doppler  period  T via 


2 sin  (?) 


It  is  important  to  note  that  there  is  statistical  bias  in  the  time-average  mean  velocity  as 
calculated  from  the  individual  velocity  data  point  realizations  from  the  output  of  the  data 
processor.  On  the  average,  the  scattering  particles  ma>  be  assumed  to  distribute  uniformly 
throughout  the  How.  Over  a long  interval  of  lime,  particles  moving  faster  than  the  time- 
average  velocity  will  be  detected  more  frequently  than  slower  velocity  particles  since  the 
probability  of  detecting  a particle  is  proportional  to  the  volume  of  fluid  swept  through  the 
probe  volume.  Reischman  and  Tiederman41  have  shown  that  an  improved  estimate  of  the 
time-average  velocity  can  be  made  by  weighting  each  velocity  realization  uxi  with  a function 
which  is  inversely  proportional  to  the  volume  flow  through  the  probe  volume  at  that  inst.  t. 
To  account  for  this  fact,  an  improved  estimate  of  the  time-average  mean  velocity  proposed  by 
Reischman  and  Tiederman41  has  been  adopted  here: 


Ux  = ux=~ 


Zk 


Here  uxi  is  inversely  proportional  to  Tj  which  is  a direct  output  of  the  DDP.  A complete 
accuracy  analysis  of  the  present  laser  Doppler  velocimeter  system  has  not  been  made. 
However,  during  many  repeated  runs  in  the  course  of  the  experiment,  it  was  found  that  the 
standard  deviations  of  the  measured  ti  were  less  than  0.02  Uu. 


' Kehdiinat  V'  M and  W <..  Tiederman.  “Lasei-Dopplei  Auetnornelei  Measurements  m Drag  Reducing 
Channel  I lows,”  Journal  of  Fluid  Mechanics,  Vol.  70,  lb  2,  pp.  369-382  (1975). 
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